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Abstract 

This work is concerned with various aspects of the formulation of the quantum in- 
verse scattering method for the one-dimensional Hubbard model. We first establish the 
essential tools to solve the eigenvalue problem for the transfer matrix of the classical "cov- 
ering" Hubbard model within the algebraic Bethe Ansatz framework. The fundamental 
commutation rules exhibit a hidden 6-vertex symmetry which plays a crucial role in the 
whole algebraic construction. Next we apply this formalism to study the SU{2) highest 
weights properties of the eigenvectors and the solution of a related coupled spin model 
with twisted boundary conditions. The machinery developed in this paper is applicable 
to many other models, and as an example we present the algebraic solution of the Bariev 
XY coupled model. 



1 Introduction 



The discovery of the quantum version of the inverse scattering method in the late seventies 
was undoubtedly a remarkable contribution to the development of the field of exactly solv- 
able models in (1 + 1) dimensions |l|]. This method provides a means for integrating models 
in two-dimensional classical statistical mechanics and (1 + 1) quantum field theory, unifying 
major achievements such as the transfer matrix ideas, the Bethe Ansatz and the Yang-Baxter 
equation. Nowadays detailed reviews on this subject are available in the literature, for instance 
seerefs. 0, |, |, i- 

We shall start this paper illustrating the essential features of this method in the context 
of lattice models of statistical mechanics. For example, consider a vertex model on the square 
lattice and suppose that its row-to-row transfer matrix can be constructed from an elemen- 
tary local vertex operator £_4j(A). This operator, known as the Lax operator, contains all 
information about the structure of the Boltzmann weights which are parametrized through the 
spectral parameter A. The operator £_4j(A) is frequently viewed as a matrix on the auxiliary 
space A, corresponding in the vertex model to the space of states of the horizontal degrees of 

L 

freedom. Its matrix elements are operators on the Hilbert space JJ^ ®Vi, where Vi corresponds 

1=1 

to the space of vertical degrees of freedom and i denotes the sites of a one-dimensional lattice 
of size L. In this paper we shall consider the situation in which the auxiliary space A and the 
quantum space Vi are equivalent. A sufficient condition for integrability of ultralocal models, 
i.e. those in which the matrices elements of the operator £^j(A) commute for distinct values 
of index i, is the existence of an invertible matrix i?(A,/i) satisfying the following property 

R{X, fi)CAiW ® CMi^i) = ® CmWR{X, yu) (1) 

where the tensor product is taken only with respect to the auxiliary space A. The matrix 
R{X, fi) is defined on the tensor product A ^ A and its matrix elements are c-numbers. A 
ordered product of Lax operators gives rise to the monodromy operator T(A) 

T(A) = CM{\)CAL-ii\) ■ ■ ■ CAiiX) (2) 
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It is possible to extend property (1) to the monodromy matrix, and such global intertwining 
relation reads 

R{\, fi)T{X) ® r(/i) = T(/i) ® T(A)i?(A, /i) (3) 

The transfer matrix of the vertex model, for periodic boundary conditions, can be written 
as the trace of the monodromy matrix on the auxiliary space A 

T(A) = Tr^nX) (4) 

From the above definition and property (3) we can derive that the transfer matrix is the 
generating function of the conserved currents. Indeed, taking the trace of equation (3) on the 
tensor A space and using the trace cyclic property we find 

[T(A),T(^)] = (5) 

Consequently, the expansion of the transfer matrix in the spectral parameter yields an 
infinite number of conserved charges. We recall that local charges are in general obtained as 
logarithm derivatives of T(A) P, 0]- Furthermore, the compatibility condition of ordering three 
Lax operators i3^i(Ai), €^^2(^2) and Cj^-^i^s) through the intertwining relation (1) implies the 
famous Yang-Baxter equation 

-R23(-^l5 A2)-Rl2(Al, A3)i?23(-^2, A3) = -Rl2(A2, A3)-R23(-^l5 A3)i?i2(Ai, A2) (6) 

where Rab{X, fi) denotes the action of matrix R{X, fi) on the spaces Va (8> H. 

Equation (3) is the starting point of solving two-dimensional classical statistical models 
by an exact operator formalism. This equation contains all possible commutation relation 
between the matrix elements of the monodromy operator T(A). The diagonal terms of T(A) 
define the transfer matrix eigenvalue problem and the off-diagonal ones play the role of creation 
and annihilation fields. The eigenvectors of the transfer matrix are constructed by applying 
the creation operators on a previously chosen reference state, providing us with an elegant 
formulation of the Bethe states. For this reason this framework is often denominated in the 
literature as the algebraic Bethe Ansatz approach. This situation resembles much the matrix 
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formulation of (0 + 1) quantum mechanics. It is well known that the harmonic oscillator can 
either be solved by the Schrodinger formalism or by the Heisenberg algebra of creation and 
annihilation operators. The later approach, however, is conceptually much simpler provided the 
relevant dynamical symmetry has been identified for a given quantum system. One successful 
example is the solution of the hydrogen atom through the 5*0(4) algebra 0. 

In this paper we are primarily interested in applying the quantum inverse scattering method 
for the one-dimensional Hubbard model. We recall that, after the Heisenberg model, the 
second one-dimensional lattice paradigm in the theory of magnetism solved by Bethe Ansatz 
method was the Hubbard model. The solution was found by Lieb and Wu in 1968 |Q using the 
extension of the coordinate Bethe Ansatz to the problem of fermions interacting via 5-functions 
[ p!0|] . Considering the success of the solution of the Heisenberg model by the inverse method 0, 
the next natural target for this program would then be the Hubbard model. However, it turns 
out that the solution of this problem followed a more arduous path than one could imagine 
from the very beginning. Indeed, nearly 18 years were to pass before it was found the classical 
statistical vertex model whose transfer matrix generates the conserved charges commuting with 



the Hubbard Hamiltonian. This remarkable step was done by Shastry |TT|, |T2|, ^ who also 
found the i?-matrix solution and thus proved the integrability of the Hubbard model from the 
quantum inverse method point of view. Shastry himself attempted to complete the inverse 
scattering program, but he was only able to conjecture the eigenvalues of the transfer matrix 
guided by a phenomenological approach which goes by the name of analytical Bethe Ansatz ||13|| . 
Subsequently Bariev presented a coordinate Bethe Ansatz solution for the classical Shastry's 
model, however on the basis of the diagonal-to-diagonal transfer matrix method [l^ . 

One of the main results of this paper is the solution of the one- dimensional Hubbard model 
by a first principle method, namely via the algebraic Bethe Ansatz approach Q For this 
purpose we will use Shastry's i?-matrix as well as the modifications introduced by Wadati and 
co-workers ||15[ . Apart from the fact that the solution of the one-dimensional Hubbard model by 
the algebraic Bethe Ansatz framework remains an important unsolved theoretical challenge in 



brief summary of some of our results has appeared in ref. flf 
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the field of integrable models, there are also other motivations to pursue this program. Recent 
developments of new powerful methods to deal with finite temperature properties of integrable 



models (see for e.g. refs. ^% ^3, |19|) show clearly that the central object to be diagonalized is 



the quantum transfer matrix rather the underlying one-dimensional Hamiltonian. The transfer 
matrix eigenvalues provide us with the spectrum of all conserved charges, a fact which could 



be helpful in the study of transport properties and level statistics behaviour |21|]. Lastly, 
there is a hope that this program is the first step towards the formulation of a general approach 
for computing lattice correlation functions 0. 

We would like to remark that the ideas developed in this paper transcend the solution of the 
one-dimensional Hubbard model. In fact, the original basis of our approach might be traced 
back to the solution of the supersymmetric spl{2\l) vertex model |2^. Very recently, we have 
shown that this method provide us with a unified way of solving a wider class of integrable 
models based on the braid monoid algebra Here we also will see that the lattice analog 



of the coupled XY Bariev chain can be solved by this technique. The unusual feature of 
the Hubbard and Bariev models is that they both have a non-additive R-matrix solution. 

We have organized this paper as follows. To make our presentation self-contained, in next 
section we briefly review the basic properties of the embedding of the one- dimensional Hubbard 
model into a classical vertex model, originally due to Shastry [|TT], [T^, |TB[. In section 3 we discuss 
the commutation rules coming from the Yang-Baxter algebra. In particular, a hidden symmetry 
of 6-vertex type, which is crucial for integrability, is found. We use these properties in section 4 
in order to construct the eigenvectors and the eigenvalues of the transfer matrix of the classical 
statistical model. The Lieb's and Wu's results as well as the spectrum of higher conserved 
charges can be obtained from our expression for the transfer matrix eigenvalues. In section 5 
we present complementary results such as extra comments on systems with twisted boundary 
conditions and a discussion on the SU (2) highest weights properties of the eigenvectors. Section 
6 is dedicated to the solution of the classical analog of the coupled XY Bariev model. Our 
conclusions are presented in section 7. Finally, five appendices summarize Boltzmann weights, 
extra commutation rules and technical details we omitted in the main text. 
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2 The classical covering Hubbard model 



We begin this section reviewing the work of Shastry [|ri|, 0, on the identification of the clas- 
sical statistical model whose row-to-row transfer matrix commutes with the one- dimensional 
Hubbard Hamiltonian. Originally, Shastry looked at this problem considering the coupled spin 
version of the Hubbard model, since in one-dimension fermions and spin-| Pauli operators 
are related to each other via Jordan- Wigner transformation. In the context of statistical me- 
chanics, however, the later representation is sometimes more appealing. Here we will consider 
the coupled spin model introduced by Shastry with general twisted boundary conditions. Its 
Hamiltonian is 

L-l 



i=l 
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where {af,af} and {t^-,t^} are two commuting sets of Pauh matrices acting on the site i of 
a lattice of size L. The second term in (7) stands for the boundary conditions cr^+i = e^^'^^af , 
rj^j^i = e^'"^'^Ti , crl^i = erf, and tI_^_i = rf where 0i and 02 are arbitrary angles < 0i, 02 < 2tt. 
The coupling constant U represents the Hubbard on-site Coulomb interaction. 

In order to relate the coupled spin model to the Hubbard model we have to perform the 



following Jordan- Wigner transformation |11 



i-l L i-l 

''fc 'i 



= n ^k<^i ' = n n ^fc^^ 



k=l k=l k=l 

where Qo- are canonical Fermi operators of spins a =t, i on site i, with anti-commutation 
relations given by {c|^, Cja'} = Sij^^y. Defining the number operator Uia = cj^Cia for electrons 
with spin a on site i and performing transformation (8) we find that 

1=1 (7=T,i i=l ^ ^ 



-e 
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where the angles 0| and 0j are given by 



0t = 01 + TT + TtNI^ , 0^ = 02 + VT + TtN^ (10) 

L 

and A^^ is the number of holes (eigenvalues of the operator ^ Qacl^) of spin cr of a given sector 

i=l 

of the Hubbard model. Therefore, the Hubbard model with periodic boundary conditions 
(01 = 01 ) is related to the coupled spin model with dynamically (sector dependent) twisted 
boundary conditions imposed. This was the reason why we started with a more general coupled 
spin model, since the two representations are fully equivalent only for free boundary conditions. 

From the point of view of a vertex model, twisted boundary conditions correspond to 
the introduction of a seam of different Boltzmann weights along the infinite direction on the 
cylinder. In practice this is accomplished by multiplying one of the elementary vertex operator, 
£^l(A) say, by a "gauge" matrix G'_4 (see section 5). Such matrix is usually related to additional 
hidden invariances of the i?-matrix Hence, although twisted boundary conditions may 

affect eigenvalues and Bethe Ansatz equations in a significative way, the relevant features of 
the integrability still remain intact. Since this section is concerned with the later point, we can 



assume periodic boundary conditions without losing generality. As Shastry ||TT|, |T^, |13| pointed 
out, the mapping of the Hubbard model (modulo above subtlety) into a coupled spin system is 
quite illuminating in searching for a "covering" vertex model. It is known that the decoupled 
spin model {U = 0) can be derived in terms of a pair of uncoupled free-fermion 6- vertex models. 
This suggests that, for the interacting model, we have to look for a copy of two free-fermion 



6-vertex models coupled in an appropriate way. Shastry |T^, determined the nature of 



this coupling by demanding that it should reproduce the higher conserved charges ||TT| |^ when 
the corresponding transfer matrix T(A) was expanded in powers of the spectral parameter A. 
The solution found by Shastry for the Lax operator is given by [|12|, |l^ 



^^.(A) = /:^,(A)/::,,(A)e'^W'^5.^i«^« (11) 



^For further discussion on Hubbard's conserved charges see refs. 
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The form of operators >C^j(A) and >C^j(A) obey the 6- vertex structure 

a(A) + 6(A) a(A)-6(A) 



and 



a{\) + h{\) a{X)-h{X) 



A 'i 



(12) 



(13) 



2 2 

where the weights a(A) and 6(A) satisfy the free-fermion condition a^(A) + 6^(A) = 1. Further- 
more, the constraint h{X) is determined in terms of the weights and the couphng U by 

U 



sinh[2/i(A)] = -a(A)6(A) 



(14) 



A second important result due to Shastry [|T^, was the solution of the Yang-Baxter 
algebra for the Lax operator (11), and thus determinating the form of the i?-matrix. The 
matrix -R(A, fi) is a 16 x 16 matrix whose non-null elements are given in terms of 10 distinct 



Boltzmann weights ai(A,/i), i 
matrix form 



10. For practical calculations it is helpful to display its 
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(15) 

where the expressions for the weights ai{X, fi) in terms of the free-fermion weights a(A), 6(A) 
and the constraint h{X) can be found in appendix A. The striking feature of this solution 
is that i?-matrix (15) is non-additive with respect the spectral parameters. In fact, after an 
unitary transformation, R{X,fi) can be written in a more compact form [0 which shows that 
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it depends on both the difference and the sum of the spectral parameters. As far we know, 
it is still an open question whether or not there exists an embedding for the Hubbard model 
satisfying the standard difference property. As a final remark we mention that an analytical 
proof that -R(A, fi) indeed satisfies the Yang-Baxter equation (6) has been recently presented 
in ref . 



We close this section presenting the graded Yang-Baxter formalism |]28[ for the Hubbard 
model. This interesting approach was pursued by Wadati and co-workers |jl5| and it has the 
advantage of making real distinction between bosonic and fermionic degrees of freedom. In the 
Hubbard model, the empty and doubly occupied sites play the role of bosonic states while the 
spin up and down states are the fermionic ones. This formalism is an elegant mathematical 
procedure^ of avoiding the subtlety on boundary condition raised in the beginning of this 
section. In other words, the graded version of the inverse scattering method guarantees that 
the "non-local" anticommutation rules of fermionic degrees of freedom is satisfied for any lattice 
sites. In general, the basic changes we need to perform is to consider the analogs of the trace 
and the tensor product properties on the graded space. For example, the graded Yang-Baxter 
for the monodromy matrix now reads 



RgiX, fi)T{\) ® Tifi) = Tifi) T{\)Rg{\ /i) (16) 

where the symbol ® stands for the supertensor product [A ® B)"^ = (— 1)^*^^-"^*^"^+^*^^^] Aac-Bbd. 
The index p{a) is the Grassmann parity of the a-th degree of freedom, assuming values p{a) = 
for bosonic specie and p{a) = 1 for fermionic ones. Other important change is on the transfer 
matrix definition, which is now given in terms of the supertrace of the monodromy matrix 

T(A) = Str^TiX) = Y.{-ir^^Ka{X) (17) 

aeA. 

There is no extra effort to obtain the matrix Rg{X, fi) from the original solution found by 
Shastry. One just have to perform a Jordan- Wigner transformation on the Lax operator (11), 
taking into account the gradation of the space of states It turns out that the graded 

•^This scheme accommodates a particular class of models having "nonultralocal" Yang-Baxter relations. For 

more general implications of nonultralocality see the recent review p9|. 
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i?-matrix is related to Shastry's solution (15) by a unitary transformation, and its explicit 
form is given by [jl5l 
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(18) 

where here we assumed that the first and the fourth degrees of freedom are bosonic = 
p(4) = 0) while the remaining ones are fermionic {p{2) = p{3) = 1). 

In the next sections we are going to use the graded formalism in order to find the ap- 
propriate commutation rules, the eigenvectors and the eigenvalues of the transfer matrix (17). 
Afterwards, we will get back to the standard quantum inverse formalism, but now with twisted 
boundary conditions. 



3 The fundamental commutation rules 

In addition to the Lax operator and the i?-matrix the existence of a local reference state is 
another important object in the quantum inverse scattering program. This is a vector \0)- such 
that the result of the action of the Lax operator on it is a matrix having a triangular form. We 
choose |0)^ = (^0) ® (oj as the standard spin up "ferromagnetic state", which in the fermionic 
language corresponds to the doubly occupied state. The action of the vertex operator in this 
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state satisfies the following property 

/^i(A)|0), 




V 








^2(A) |0), 









(19) 



^2(A) |0), t 

^3(A)|0)J 

where the symbol \ represents arbitrary non-null values and the functions a;i(A), a;2(A) and 
(A) are given by 



,;i(A) = [a(A)]VW , o;2(A) = a(A)6(A)e-^W , ^3(A) = [6(A)]V(^) 



(20) 



The global reference state |0) is then defined by the tensor product |0) = J]^® \0)-. This 

i=l 

state is an eigenstate of the transfer matrix since the triangular property is easily extended to 
the monodromy matrix. In order to construct other eigenstates it is necessary to seek for an 
appropriate representation of the monodromy matrix. By this we mean a structure which is 
able to distinguish creation and annihilation fields as well as possible hidden symmetries. The 
triangular property of the Lax operator suggests us the following form 

/B{X) B{X) F{X) \ 



r(A) 



C(A) i(A) B*{\) 
\C{X) C*{X) D{X)J 



(21) 



4x4 



where B{X), C*{X) and B*{X), C{X) are two component vectors with dimensions 1x2 and 
2x1, respectively. The operator A{X) is a 2 x 2 matrix and we shall denote its elements by 
Aab{X). The remaining operators -B(A), C(A), D{X) and F{X) are scalars. In this paper we 
will use the symbol ABCDF to refer to the above way of representing the elements of the 
monodromy matrix. We recall that such Ansatz is quite distinct from the traditional ABCD 
form proposed originally by Faddeev and co-workers 0, 0, 0]. 

In the ABCDF representation the eigenvalue problem for the graded transfer matrix be- 
comes 

[5(A) - Aaa{\) + D{X)] 1$) = A(A) 1$) (22) 



a=l 
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where A (A) and |$) correspond to the eigenvalues and to the eigenvectors, respectively. As a 
consequence of the triangular property we can derive important relations for the monodromy 
matrix elements. For the diagonal part of T{X) we have 

5(A) |0) = [a;i(A)]^ |0) , D(A) |0) = [^3(A)]^ |0) , ia,(A) |0) = [a;2(A)]^ |0) for a = 1, 2 (23) 

— * — # 

Also one expects that the operators -B(A), B*(\) and F(X) play the role of creation fields 
over the reference state |0). It also follows from the triangular property the annihilation 
properties 

C(A)|0) = , C*{X)\0) = , C(A)|0) = , iafe(A) |0) = for a 7^ b (24) 

To make further progress we have to recast the graded Yang-Baxter algebra in the form 
of commutation relations for the creation and annihilation fields. In general it is not known 
how and when such job can be performed for a particular representation, and one could surely 
say that the "artistic" part of the algebraic Bethe Ansatz construction begins here. Within 
the ABCDF formalism, the solution of this problem turns out to be more complicated than a 
similar situation occurring for the 6-vertex model [|I], 0, ^ and its multi-state generalizations 
31]. The new feature present here is that we have a mixture of two classes of creation 



fields, the non-commutative vectors B{X) or B*{X) and one commutative operator represented 
by F{X). We shall start our discussion by the commutation rule between the fields -B(A) 
and B{fi). In this case the relation that comes out from the Yang-Baxter algebra is not the 
convenient one for further computations. It turns out to be necessary to perform a second step 
which consists in substituting the exchange relation for the scalar operators B{X) and F{fi) 
(see equation (38)) back on the original commutation rule we just derived for the fields B{X) 
and B{fi). The basic trick is to keep the diagonal operator -B(A) always in the right-hand side 
position in the commutation rule After performing this two step procedure we are able 
to get the appropriate commutation rule, which is 

5(A) ® S(/x) = ^^^Mfi) ® S(A)].f(A, /i) - ^^^j^{F{X)Bi^^) - F(/x)5(A)K (25) 



11 



— * 

where ^ is a 1 x 4 vector and f(A, //) is an auxiliary 4x4 matrix given by 

/I 0\ 

a{\, fi) b{X, fj,) 
6(A,/i) a{\, fi) 

Vo ly 

and the functions a(A, /i) and 6(A, //) are given in terms of the Boltzmann weights by 



e= (0 1 -1 0) , r(A,//) 



(26) 



a{X,ii) 



Q!3(A, ii)a7{X, fi) + alo{X, n) 



6(A,//) = 



q;6(A, iJ^)a7{X, fi) + q;?o(A, A*) 



(27) 



ai{X, fi)a7{X, fj,) ai{X, fi)a7{X, n) 

It turns out that the auxihary matrix f(A,/i) is precisely the rational i?-matrix of the 
isotropic 6-vertex model or the XXX spin chain. In order to see that, we first simplify a bit 
more the auxiliary weights a(A, //) and 6(A, //) with the help of identities (A.10-A.12). We find 
that they satisfy the following relations 

as{X, n)ag{X, n) 



a(A, //) = 1 - 6(A, //) , 6(A,/i) = 



q;i(A, /i)Q;7(A, /i) 



(28) 



Next we simplify as much as possible the ratios and in terms of the free- 

fermion Boltzmann weights and the constraint h{X). After some algebra we write these ratios 



as 



ag{X,fi) _ a(/i)6(A)e^[^(^)-^Wi - a(A)6(/x) 
ai{X,fi) ~ 6(A)6(At) + a(A)a(/x)e2[MM)-MA)] 

as{X,i2) _ b{X)b{i2) + a(A)a(^)e^[^(^)+^(^)] 
a7{X,fi) ~ a(//)6(A)e2[M/^)+MA)] - a{X)b{fi) 

Now if we take into account the identity 



a(A) 



,2hi\) 



6(A) 



-2h(X) 



a(A) 



-2h{X) 



6(A) a(A) 6(A) 

and perform the following reparametrization 



a(A)' 



,2hiX) 



U 



X 



a(A) 



,2h{X) _ 



6(A) 



2/i(A) _ 



U 



6(A) a(A) 
we finally can rewrite the auxiliary weights as 

U 



a(A,/x) 



jl-X + U 



, 6(A,/i) 



/i — A 
/i- A + C/ 



(29) 
(30) 

(31) 
(32) 
(33) 
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Clearly, these are the non-trivial Boltzmann weights of the isotropic 6-vertex model. This 
is an important hidden symmetry, which is known to play a decisive role on the exact solution 
of the Hubbard model since the work of Lieb and Wu ||^. The derivation of this symmetry in 
the context of the quantum inverse scattering program is however a rather non-trivial result. 
One of the virtues of this result is that it becomes valid for the generator of the commuting 
conserved charges and not only for the Hubbard Hamiltonian. Moreover, we also recall that 
this symmetry is of relevance to the Yangian invariance of the Hubbard model which emerges 
in the thermodynamic limit ||3^, |33[| . 

To solve the eigenvalue problem (22) we still need the help of several other commutation 
relations. For instance, the commutation rules between the diagonal and creation operators 
play an important role in the eigenvalue construction. It turns out that in some cases we have 
to take into account similar trick discussed above. This is specially important for the field 
A{X), where we have to use an auxiliary exchange relation between the operator B{fi) and 
B*{X), in order to obtain a more appropriate commutation rule with the creation operator 
B{X). In general, the task is quite cumbersome and here we limit ourselves to list the final 
results. The commutation relations between the diagonal fields and the creation operator B{X) 
are 

i(A)®5(/i) = -z ^^S^'^l [%)^i(A)].f(A,/i) + ^ "'S^^| ^(A)0i(/i) 

ag{X,n) a9{X,fi) 



^ aio(A,/i) 
a7(A,/i) 



ag[X,fl) ag[X,fi) 



(34) 

i?(A)5(/.) = ^^^B{^^)BiX) - ^^^B{X)Bi^,) (35) 

D{X)BM = -,^^B(,)D(X) + ^^F(u)C'(X) 

"^^^'^^F(A)C*(/i) - t^4^l[B*iX) ® Ai^^)] (36) 



a'i{X, jj,) a'j{X^ jj) 

while those for the scalar field F{X) are 

i.,(A)F(/.) = [1+ ^f^^'^^ ]F(/.)i,,(A) - , F{X)AM 
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+^^4^^[^(A) B*(^)U + i^^AB*(X) (37) 
ag{X,iJ,) a;8(A,/x) 

B(A)F(^) = ^^F^^)B(X) - '^F(X)BM + ® BM}-? (38) 

a7[ii,X) a7[ii,X) a7{ii,X) 

DWF(,) = ^^F(^)D(X) - ^^F(X)D(^) - ^^^H^^f {ir-(A) » B-(^)} (39) 
q;7(A,//) Q;7(A,//j Oi7{X,ii) 

-*, -* 
where ^ stands for the transpose of ^. Furthermore, the relations closing the commutation 

rules between the creation operators B[X) and F[X) are 

[F(X),F(p)] = (40) 

Finally, it remains to consider the commutation rules for the creation field B*{X). To avoid 
overcrowding this section with more heavier formulae we have collected them in appendix 

— * 

B. We see that they are quite similar to those we just derived for the field B{X). In fact, 
it is possible to establish an equivalence between these two sets of commutation rules if we 
formally interpret the symbol * as a mathematical operation acting on the elements of the 
monodromy matrix. For lack of a better name we call it "dual" transformation and we impose 
that it satisfies the following properties: (0*(A))* = 0(A), A*{X) = -A\X), B*{X) = D{X), 
F*{X) — F{X) and C*(A) = C{X). Applying the "dual" transformation on the commutation 
rules of field -B(A) we obtain those for the field B*{X) with new Boltzmann weights a* (A, /x; h) = 
aj{X, fi; —h), where, for sake of clarity, we stressed the dependence on the constraint h{X). 
This means that the functional form of the weights remains unchanged but now we have to 
perform the transformation h{X) — > —h{X) {U — > —U). We recall that in this last step we 
used the following identities for the Boltzmann weights: aj{X, h) — CKj{iJ,, A; —h) j — 1, . . .7, 
as{X,fi;h) = ag{fi, X; —h) and aio{X,fi;h) = — q;io(/x. A; — /i). Therefore, we expect that the 
construction of the eigenvectors will be based either on the pair of fields B[X) and F{X) or on 
the "dual" ones B*{X) and F{X) rather than on a general combination of the three creation 
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fields. This redundance is in accordance to what one would expect from the space of states of 
the Hubbard Hamiltonian, since at a given site we can either create a single electron (spin up 
or down) or a pair of electrons with opposite spins. We remark that such "duality" property is 
not particular to the Hubbard model but it is rather a general feature present in our framework 



At this point we have set up the basic tools to start the construction of the eigenvectors 
of the eigenvalue problem (22). In next section we will show how this problem can be solved 
with the help of the commutations rules (25 — 26; 34 — 42) and few other relations presented 
in appendix B. 

4 The eigenvectors and the eigenvalue construction 

The purpose of this section is to solve the eigenvalue problem for the graded transfer matrix. 
We shall begin by considering the construction of an Ansatz for the corresponding eigenvectors. 
The multi-particle state are going to satisfy an important recurrence relation. We will see that 
the eigenvalue problem (22) has a nested structure, i.e. it will depend on the solution of an 
inhomogeneous auxiliary problem related to the 6- vertex hidden symmetry. 

4.1 The eigenvalue problem 

The eigenvectors of the transfer matrix are in principle built up in terms of a linear combination 
of products of the many creation fields acting on the reference state. These Bethe states are 
often thought as multi-particle states, characterized by a set of rapidities parametrizing the 
creation fields. Before embarking on the technicalities of the construction of an arbitrary 
n-particle state we first define it by the following scalar product 



m 




(43) 



where the mathematical structure of vector 




creation fields. At this stage the components of vector JF are simply thought as coefficients of 
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an arbitrary linear combination which would be determined later on. This reflects the "spin" 
degrees of freedom of the space of states and we shall denote such coefficients by ^""•■■"i where 
the index at run over two possible values = 1,2. 

Let us now turn our attention to the construction of vector $„(Ai, . . . , A„). As mentioned at 
the end of the previous section, it is sufficient to look for combinations between the fields B{X) 
and F{X). In general, there is no known recipe which is able to provide us with an educated 
Ansatz for this vector and as it is customary we shall start the construction considering few 
particle excitations over the reference state. A single particle excitation is made by creating 
a hole of spin up or down on the full band pseudovacuum |0). From the point of view of the 
inverse scattering method this excitation is represented by ^i(Ai) = -B(Ai) and consequently 
the one-particle state is 

|$i(Ai)) = B{X,).f |0) - 5a(Ai).F« |0) (44) 

where from now on we assume sum over repeated index. 

It is not difficult to solve the eigenvalue problem (22) for such one-particle state. If we use 
the commutation relations (34-36), and the pseudovacuum properties (23-24) we find that the 
one-particle state satisfies the following relations 

B{X)\M>^^))^z^^^^M{Xr\M^,))-t^^^ (45) 
tt9(Ai,A) a9(Ai,A) 

a8(A, Ai) ^ (1) \ -"lol-^^^Or, , ^^ Mir 



D(X) |$i(AO) = -^^=^^hnrT[^3(A)]^ ^i(Ar) - ^ ; MXirUB*(X) I)].:F |0) 

Q!7(A,Ai) / q;7(A,Ai) 



(46) 



EAaa(A)|$i(AO) = 

a=l "91^5 ^ij 



ag{X, Aij 

- Cixx'^ M>^immx) ^ |Q) (47) 

"tIA, Aij 

where / is the 2x2 identity matrix. The terms proportional to the eigenvector |$i(Ai)) are 
denominated wanted terms because they contribute directly to the eigenvalue. The remaining 
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ones are called unwanted terms and they can be eliminated by imposing further restriction on 
the rapidity Ai. This constraint, known as the Bethe Ansatz equation, is given by 

_a;2(Ai)_ 

It is now straightforward to go ahead and to determine the one-particle eigenvalue. However, 



it is convenient to start introducing suitable notation which can be extended to accommodate 
multi-particle states. With this in mind, we define the following auxiliary eigenvalue problem 

r(i) (A, = (A, Ai)^"^ = A(^) (A, Ai)^^i (49) 

and we see that, in terms of equation (49), the one-particle eigenvalue can be expressed by 

a9(Ai,Aj q;7(A, Aij q;9(A, Aij 

Up to the level of the one-particle state there is no extra effort to solve the corresponding 
auxiliary problem. Considering the 6-vertex structure of matrix f(A, //) it is easily seen that 
the solution is 

A(i)(A,Ai) = l + 6(A,Ai) (51) 

We next turn to the analysis of the two-particle state. We expect that such state will be a 
composition between two single hole excitations of arbitrary spins and a local hole pair with 
opposite spins. The former is made by tensoring two fields of B{\) type while the later should 
be represented by F{\). The vector ^ has also a physical meaning. It plays the role of an 
"exclusion" principle, forbidding two spin up or two spin down at the same site. Thus, an 
educate Ansatz for the two-particle vector should be the linear combination 

^2(Ai, A2) = S(Ai) (8) B{\2) + |F(Ai)S(A2)#(Ai, A2) (52) 

where (Ai, A2) is an arbitrary function to be determined. We found also convenient to add 
the diagonal field -B(A2) on the right-hand side of the two-particle vector Ansatz. We see that 
when the Ansatz (52) is projected out on the subspace of equal spins, no contribution coming 
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from F{X) appears, which is in perfect accordance to what one would expect from the Pauh 
principle. In other words, using the definition (43) we have 

|$2(Ai, A2)) = Bi{X,)Bj{X2)J^^'\0) + [u;i(A2)]^F(Ai)#(Ai, A2)(^'' -^'') (53) 

In order to tackle the eigenvalue problem for the two-particle state, besides the commutation 
rules of the last section, we have to use extra relations between the fields -B(A), -B*(A), C{X) 
and C*{X). These relations have been summarized in the beginning of appendix B. After 
turning the diagonal fields over the two-particle state, we find that there are two classes of 
unwanted terms. The first class we call "easy" unwanted terms because they are only produced 
by the same diagonal operator ( A{X) or D{X) ) and they can be eliminated by an appropriate 
choice of function ^o^^(Ai, A2). There are three terms of this sort 

F{X)D{X,)B{X2) , B{X).B*{X^)B{X2) , l[B*{X) ^ B*{X^)]B{X2) (54) 

and all of them are cancelled out provided we chose function ^o^'*(Ai, A2) as 

9f(Ai,A,) = .^^HS^ (56) 
a7[Xi, X2) 

Now, besides the wanted terms, we arc only left with standard unwanted terms, i.e. those 
that require further restriction on the rapidities. We shall see below that these terms can be 
simplified in rather closed forms with the help of the two-particle auxiliary problem. Similar 

to the one-particle analysis, the auxiliary eigenvalue problem is figured out by looking at the 

2 

wanted terms coming from the operator ^ Aaa{X). Considering the commutation rule (34) we 

a=l 

soon realize that the two-particle auxiliary problem is 

With the above information we move on simplifying as much as possible the action of the 
diagonal fields on the two-particle state. We keep in mind that we want to present the results 
in a way that would be amenable to multi-particle states generalization. After a cumbersome 
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algebra we find that 



S(A)|$2(Ai,A2)) 



-E[^i(A,)]"|*S'^(A,A,;{AJ)) 

+H,{\ Ai, A2)[u;i(AiVi(A2)]^ ^i'^A, A,-, A;; {A^}) 



(57) 



L'(A)|$2(Ai,A2)) = 



MA)]^n-z^#^l*2(Ai,A2)) 



"7(A, Aj) 



Y^[u;,iXjrA('\X = A„ {A,}) \¥^\X, A,; {A,})) 



+i/2(A, Ai, A2)[^2(Ai)w2(A2)]^ (A, A,-, A^; {AJ) 



(58) 



^i„„(A)|$2(Ai,A2)) 



^ .ai(A, Aj^ 



A«(A,{A,})|$2(Ai,A2)) 



0=1 



[-^^^^n-^.,A,A.) 

- X:[a;2(A,)]^A«(A = A,, {A,}) |*«(A, A,; {A,})) 

-E['^i(A,)]"|^S'^(A,A,;{A,})) 

+Hs{X, Ai, A2)[a;i(Ai)u;2(A2)]^ |*?Ha, A,-, A^; {AJ)) 

+i/4(A, Ai, A2)[a;i(A2)a;2(Ai)]^ \^^^\x, X„ Xf, {AJ)) (59) 

For sake of clarity we have shortened the notation for the unwanted terms and represented 
them by the eigenfunctions '^i\X, Xj; {Xk})^ and ^o^^(A, A^; {A^D^ We see that there are 
three classes of unwanted terms and their explicit expressions in terms of the creation fields 



are 



^?{X,X,-{Xi}) 



,. «5(Aj, A) -p|- . Q2(Afc, Xj) 

<^9{Xj,X) ag{Xk,Xj] 

k^3 



5(A) ® $i(Afc)0«(A,-; {A,}).^ |0) (60) 



,(2)m \ r\.l^\ _ „-"io(A, Aj; -j-j- ^ 



M'^(A,A,;{A,}) 



«2(A,, A,) ^^^^^^^^^ ^ ^^^j ^ $i(A,)6f (A,; {A,})./ 10) 



«7(A,Aj) Q;9(Afe,Aj) 



(61) 
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Before proceeding with a discussion of the resuhs, we should pause to comment on the 
"brute-force" analysis we performed so far for the two-particle state problem. Roughly speak- 
ing, one can estimate the wanted terms by keeping the first term of the right-hand side of the 



This procedure gives us the coefficients proportional to the first part of the eigenvector and 
to show that this is also true for the second part wc need to use some identities between the 
Boltzmann weights. The situation for the unwanted terms is even worse due to the proliferation 
of many different terms, common in a Bethe Ansatz "brute-force" analysis. The "ordering" 
factor just accounts for these many different contributions to the unwanted terms. Later on it 
will become clear that the origin of this factor is due to a permutation property satisfied by the 
two-particle eigenvector. In appendix C wc provide the details about the less straightforward 
simplifications carried out for the two-particle state, since some of them will be also useful 
to multi-particle states as well. Finally, within a "brute force" computation, we have found 
nine contributions to the third unwanted term which come from many different sources. It is 
possible to recast them in terms of four functions Hi{x,y, z) i = 1, . . . ,4: whose expressions are 



commutation rules (34-36) when we turn the diagonal fields over the creation operators B{Xj). 



H2{x,y,z) 



H4{x,y,z) 



Hi{x,y,z) 



H3{x,y,z) 



a2{y,x)a5{z,x)aiQ{y,x) _ . a4{y,x)aio{y, z) 
ag{y, x)ag{z, x)a'i{y, x) aj{y, x)a-j{y, z) 
a5{x,y)aiQ{x,z) _ . aA{x,y)aiQ{y, z) 
aj{x,y)aj{x,z) a7{x,y)a7{y, z) 
aio{x,y)a5{x,y)a5{y, z) _ ^^ a2ix,y)a5{x, z)aio{x,y) 
a7{x, y)ag{x, y)ag{y, z) ag{x, y)ag{x, z)a'!{x, y) 
. aw{x,y)a5{x,y)a5{y,z) _^ . ai{x, y)aio{x, z)a5{x, y) [1 + a{x, y)] 
a7{x, y)ag{x, y)ag{y, z) ag{x, y)a7{x, z)a8{x, y) 



al{x,y)aio{y,z) 



(64) 



a8{x,y)ag{x,y)a7{y,z) 
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Now we return to the discussion of the two-particle state results. For the first two classes of 
unwanted terms we only have two main contributions and from equations (57-59) it is direct to 
see that they vanish provided that the rapidities satisfy the following Bethe Ansatz equations 

L 

= A«(A = A,,{A,}), i = l,2 (65) 
Furthermore, the above Bethe Ansatz equations are also sufficient to cancel out altogether 

— * — * 

the four contributions proportional to the unwanted term F{X)^.T. A simple way of seeing 
that is first to factorize a common factor [uj2{Xi)uj2{X2)]^ for all the four terms. This is done 

by substituting the values [cJi(Ai)]'^ and [c<Ji(A2)]^ given by the Bethe Ansatz equations (65) 
and by using the following two-particle relations 

A^'\X = X,,{Xi})Cf = [KX,,X2)-aiX,,X2)]Cf 

AW(A = A2,{Ai})e.J? = [6(A2,Ai)-a(A2,Ai)]e.J? (66) 

After putting all these simplifications together, one is still left to verify that the following 
identity 

Hi{x,y,z) + H2{x,y,z) = H3{x,y, z)[b{y, z) - a{y, z)] 

+H4{x,y,z)[b{z,y) -a{z,y)] (67) 

is satisfied. At this point we note that there is a way of rewriting the term H4{X, Xi, X2) in 
a more symmetrical form. This technical point is discussed in appendix C and proved very 
useful in carrying out the cancellation mechanism for general multi-particle states. 

Finally, from equations (57-59) we can read directly the wanted terms, and the two-particle 
eigenvalue is 

A(A,{A.}) = ^(Arn,^4^+MArn-.^4^ 

-MX)]'{l-^^^A^'\X,{XJ}) (68) 

Now we reached a point which is typical of nested Bethe Ansatz problems, i.e. the solution 
of the two-particle auxiliary problem is no longer trivial and it is necessary to implement a 
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cc^i(Ai) 

<^2(Aj) 



second Bethe Ansatz. We will postpone this discussion until the next subsection in which 
we will present the solution of this problem for general multi-particle states. Although, for an 
integrable model, it is beheved that the two-particle sector contains the essential features about 
the general structure of the eigenvalues and the Bethc Ansatz equations, similar situation for 
the eigenvectors is still less clear. Before considering this problem, it is wise to look first for 
an alternative way of starting with a general Ansatz, since a "brute force" analysis proved 
to be rather intricate even for the two-particle state. In fact, there is a symmetry which we 
have not yet explored. It consists of seeking for eigenvectors which are in some way related 
to each other via permutation of the rapidities. This idea goes along the lines the usual 
pseudomomenta symmetrization imposed to coordinate Bethe Ansatz wave functions. For 
example, let us consider the two-particle vector in which the constraint ^q^^(Ai, A2) has been 
fixed as in equation (55). Then, it is possible to verify that the following exchange property 

^2(Ai, A2) = ^4t^*2(A2, Ai).f(Ai, A2) (69) 
q;2(Ai, A2) 

is satisfied. In order to show that, we used a remarkable relation between vector ^, the auxiliary 
matrix f(A, /i) and the Boltzmann weights given by 

^./x X aio(A,/i)a7(M, A)a2(A,^) 

C-?^A,/i = — - — ^ — - — - — -i 70 

a^[X, iJ,)aio{ii, A)Qii(A, //) 

Alternatively, we can reverse the arguments demanding that the eigenvectors satisfy the 

exchange symmetry (69). This procedure gives us a restriction to function ^q^^(Ai, A2) and it 

is an elegant way of fixing the linear combination from the very beginning. Now it is easy 

to understand the reason why an "ordering" factor had emerged in the "brute-force" analysis 

of the two-particle state. For example, the simplest way to generate the unwanted terms 

B{X) (g) $i(Ai) and [i.{B*{\) (g) /)] (g) $i(Ai) is by using the right-hand side of equation (69) 

instead of the left-hand side we used in the whole "brute force" analysis. In this way we 

obviously generate only one contribution to such unwanted terms which carries the "ordering" 

factor explicitly. 

In principle, such symmetrization mechanism can be implemented to any multi-particle 
state, and as we shall see below, it indeed help us to handle the problem of constructing a 
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general n-particle state Ansatz. We will start our discussion considering the three-particle 
state. This state is expected to be a composition between the term representing the creation 
of three holes (arbitrary spins) on different sites and the three possible ways of combining pairs 
of holes with a single excitation. Within our algebraic framework the Ansatz encoding these 
features is 

$3(Ai, A2, A3) = ^(Ai) 5(A2) ® ^(A3) + [5(Ai) ®e>(A2)5(A3)]#(Ai, A2, A3) 

F(Ai)S(A3)S(A2)]#(Ai, A2, A3) + ® F(Ai)S(A2)S(A3)]#(Ai, A2, A3) 

(71) 

where the coefficientes gj^\Xi, A2, A3) are going to be determined assuming a priori an exchange 
property (cf. equation (78)) for the Ai <-> A2 and A2 ^ A3 permutations. To see how this works 
in practice, let us first implement the permutation between the rapidities A2 and A3. To this end 
we use the commutation relation (25) to reorder these rapidities in the permuted three-particle 
vector ^3(Ai, A3, A2). This allows us to write the following relation 

^4^^$3(Ai, A3, A2).f23(A2, A3) = S(Ai) ® B{X,) B{Xs) + i^^4^^[B{X^) ® iF{X,)B{X,)] 
«2(A2,A3) q;7(A2,A3) 

I roM \^cTpr\ \ Mr •'^io(A2,A3) q;i(A2,A3) 
+ [B{Xi) ^F{Xs)B{X2)][-t — + 

Q;7(A2,A3j Q;2(A2,A3j 

X^O^^(Ai, A3, A2).f23(A2, A3)] 

+[e ® F(Ai)5(A2)5(A3)]^4r44^i'^('^i' A2).r23(A2, A3) 

q;2(A2, A3) 

+ [e ® F(A05(A3)5(A2)]^4r44^?^(^i' A2).r23(A2, A3) 

Qi2(A2, A3) 

(72) 

Imposing the exchange property to the thrcc-particlc state, i.e. that the right-hand sides of 
equations (71) and (72) are equal, we are able to derive constraints to functions Qj (Ai, A2, A3). 
We find that it is sufficient to have 

-(3)/\ A \ \ .ttlo(A2, A3) 

90 Ai, A2, A3) = ^ — , . 73) 
a7(A2, A3) 
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and 

g\X„ X„ A3) = '''i^'M #(Ai, A3, A2).f23(A2, A3) (74) 

"2IA2, Asj 

where we used the identities r23(A2, A3).f23(A3, A2) = / and ai(A, n) — Q;2(/i, A). We recall that 

relation (70) helps us to cancel out the third term of equation (72). Now it remains to determine 

function ^^^(Ai, A2, A3) and this can be done by using the permutation between the variables 

Ai and A2. The technical steps of this computation are more involving, since it is necessary to 

use other commutation rules and some identities between the Boltzmann weights. The details 

are presented in appendix D and here we quote our result for the remaining functions 

^(3)/> > , ^ .Q;io(Ai, A2) .Q;2(A3, A2) 

fifl XAi, A2, A3) = I — - — —z- 



a7(Ai,A2) a9(A3,A2) 
-(3)/x X A\ .aio(Ai, A3) .ai(A2, A3) ^ 

92 {^1,^2, h) = t — 7T — rY^^\ — r^r-23(A2, A3) (75) 

«7lAi,A3) a9(A2,A3j 

To make sure we are on the right track, we have checked that the three-particle "easy" 
unwanted terms are automatically canceled out provided we fix the constraints gf\Xi, A2, A3) 
as in equations (73) and (75). We note that functions g^\x,y,z) and g^\y-,z) are identical, 
and this allows us to rewrite the three-particle vector in terms of the following recurrence 
relation 



$3(Ai,A2,A3) = 5(Ai) ® $2(A2,A3) 

+ ^ [i® F(Ai)$i(A2, . . . , A,-i, A,+i, . . . , A3)i^(A,)] gfl,{X,, A2, A3) 



(76) 



This expression is rather illuminating, because it suggests that we can write a general n- 
particle state in terms of the (n— l)-particle and (n — 2)-particle states via a recurrence relation. 
Prom our expressions for the two-particle and the three-particle states it is not difficult to guess 

that the n-particle vector should be given by 

$„(Ai,...,A„) = S(Ai)®$„_i(A2, ...,A„) 

n 

+ E i^(Ai)4-2(A2, . . . , A,-i, A,+i, . . . , A„)5(A,)] ^(:i\(Ai, . . . , AJ 

(77) 
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where here we formally identified $o with the unity vector. Our next step is to implement the 
symmetrization scheme for such multi-particle state Ansatz. The best way to proceed here is 
to use mathematical induction, i.e we assume that the {n — 2)-particle and the {n — l)-particle 
states were already symmetrized to infer the constraints ^^"•'(Ai, . . . , A„) for the n-particle state. 
For this purpose we impose that any consecutive permutation between the rapidities Aj_i and 
Xj {j — 2, . . . ,n) satisfies the following exchange property 

^„(Ai, . . . , Xj-i, Xj, . . . , A„) — — — ^--^$„(Ai, . . . , Xj, Aj_i, . . . , A„).fj_ij(Aj_i, Xj) (78) 

a2[Aj-l, Aj) 

where the indices under fj_ij{Xj^i, Xj) emphasize the positions on the n-particle space 1 ® 
... (8) j — 1 J ... ® ^ in which this matrix acts non-trivially. 

Now starting with the latest permutation j — n we go ahead comparing the terms pro- 
portional to [^0 F(Ai)$n-2(A2, . . . , Aj_2, Xj, Xn)B(Xj-i)] and [f ® F(Ai)C-2(A2, . . . , Aj_i, 
Aj+i, . . . , Xn)B{Xj)] in both sides of the exchange relation (78). At each step, this yields a set of 
relations between the functions gj^\Xi, . . . , A„) which are further simplified by using explicitly 
both the unitarity condition and the Yang-Baxter equation for the auxiliary r-matrix. Up to 
J = 3 we find that such functions satisfy the following recurrence relation 



^j"\(Ai, . . . , Xj-i, Xj, . . . , A„) — ^' ,^}gj%{Xi, . . . , Xj, Aj_i, . . . , A„)fj_ij(Aj_i, Xj) (79) 

a2[Aj-i,Aj) 



Next we implement the symmetrization Ai ^ A2 along the lines sketched in appendix 
D for the three-particle state. In this case we have to eliminate the term proportional to 
<H) -F(Ai)$„_2(A3, . . . , A„)S(A2)] which only occurs in the left-hand side of the exchange 
relation (78). This condition helps us to determine the expression for the first constraint and 
we have 

^(n)/^ ^ N . «lo(-^l, ^2) /V . «2(Afc, A2) , . 

g\ '[Xl,...,Xn) =^ ^ X 11^ ^ /^ ^ X 80 

Q;7(Ai,A2) ^J'g Q;9(Ajfc,A2) 

Finally, the set of relations (79) and (80) are solved recursively and we find that the n- 
particle vector is 

$„(A„...,A„) = 5(AO0^„-,(A2,...,A„) + t.^^Hi^ 

j=2 (^7[M,Xj) (^9[Xk,Xj) 

25 



X 



® F(Ai)$„-2(A2, • • • , Aj_i, Aj+i, . . . , Xn)B{Xj 



X n ^Y^^rk,k+i{Xk, Xj) (81) 
At this point it is fair to remark that the recursive way we found for the eigenvectors were 



inspired to some extent on an early work of Tarasov on the Izergin-Korepin model . Our 
construction, however, has the important novelty of allowing a general "exclusion statistics" 
between the non-commutative and the commutative creation fields and therefore paving the 
way for further applications and extensions. Indeed, the non-trivial way that both the "ex- 
clusion" vector and the auxiliary matrix enters in the eigenvectors expression (81) makes our 
formula rather general, being able to accommodate the solution of a wider class of integrable 
models. This situation has to be contrasted to that of multi-state 6-vertex generalizations 



31], in which the eigenvectors are easily given by tensoring the creation fields and there is 



no explicit dependence of the underlying algebra. 

Let us now return to the problem of finding the eigenvalues of the transfer matrix, keeping 
in mind the recurrence relation (81) for the eigenvectors. To gain some insight about this 
problem we first investigate how the wanted and unwanted terms are collected for the three- 
particle state. Besides the commutation rules for the diagonal fields, we also have to use our 
previous results for the two-particle state (cf. (57-59)) wherever there is the need to carry 
the diagonal operators through the vector $2(^2,^3). This recursive way not only helps us 
to better simplify the wanted terms but also makes it possible to gather the unwanted terms 
in rather closed forms. This analysis is presented in appendix D since it still involves some 
extra technicalities. Having at hand the two-particle and the three-particle data we can move 
forward to the analysis of the four-particle state and so forth. In general, for n < 3, the 
knowledge of the {n — l)-particle and the {n — 2)-particle results dictates the behaviour of 
the n-particle state. By using mathematical induction we are able to determine the general 
structure for the multi-particle states and the final results are 

i?(A)|<|.„(Ai,...,A„)) = [^i(A)]^nz "^i^^'^| |1'„(Ai,...,A„)) 

f-Jl ag{Xj,X) 
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-^[a;i(A,)]^|vl>«,(A,A,;{A,})) 

n j — 1 

+ A,)[u;i(A,Vi(A,)]^ |xI/(;^22(A, A„ A^; {A,})) (82) 

j=2 1=1 



D(A)|$„(Ai,...,A„)) = [u;3(A)]^n 



^1 q;7(A, Aj) 



|$„(Ai,...,A„)) 



- E[^2(A,)]^A«(A = A,, {A,}) ¥:UK a,; {A,}); 
+ E E ^2(A, A;, A,)[a;2(Ai)a;2(A2)]^A(^) (A = A„ {A^}) 

3=2 1=1 

xA«(A = A^, {Afe}) *£l2(A, A,-, A^; {A,}); 



(83) 



Eiaa(A)|$n(Ai,...,A„)) - [u;2(A)]^ J] 



a=l 



: Ql(A, Aj) 
q;9(A, Ay) 



AW(A, {AJ) |*„(Ai,...,A„)) 



E[c^2(A,)]^A(i)(A = A„ {A,}) <^2i(A, A,; {A,}) 



-Eki(A,)]^ ¥:Ux,X,;{Xi})] 

- EE ^3(A, A,, A,)[a(A,, A,) - 6(A,, A,.)][u;i(A,)a;2(A,)]^ 

j=2 1=1 

xA«(A = A„ {A,}) \j, Xi; {A,})) 

- E E ^3(A, A,, Xi)^^^^MXj)u;,{X,r 

j=2l=l Oi2[Al,Al) 



xA(i) (A = A,, {A,}) A,-, A,; {A,})) 



(84) 



Similarly to what happened to the two-particle and the three-particle cases we have three 
families of unwanted terms. As before they are written in terms of the creation operators and 
the general expressions are 



n^i(A,A,;{A,})) = 



,(1) f\ \ .i\A\\ - . «5(Aj,A) " . a2(Afc,Aj) 

"9(Aj, A) l}^ ag{Xk, Xj) 

xdf\xr,{x,}).:F\o) 



5(A) (8)$„_i(Ai,...,Aj,...,A„) 



(85) 
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€U\ Ay; {A,})> = i^^^^ n '^^^^AUB'W » /)] » *„-:(A., .... A„ .... A„) 



«7(A,Aj) a9(Afc,Aj) 
AW 



xC>.^^(A,;{A4).^|0) (86) 



¥:UX, A„ A,; {A.})) = n ^^^r^T^^^^r^T^^^)^''® '^"-^(^i' • • • , A,, . . . , A„ . . . , A„) 

xdlf{Xi,Xr,{x,}).f\o) (87) 

where the symbol Xj means that the rapidity Xj is absent from the set {Ai,...,A„}. For 
n > 3 it was necessary to introduce a second "ordering" factor in order to better represent the 
third type of unwanted terms (cf. appendix D). Its task is similar to that played by the first 
"ordering" factor with the difference that now two rapidities are reordered. In other words, 
this second "ordering" factor brings the rapidities A; and Xj {I < j) to the first two positions 
in the eigenvector formula (81), and a simple calculation shows that its expression is 
/n(2)/\ \ r\ ^^ tt <^i(Afc, Aj) ^ \ \ culAfcjAj)^ , , 

^Ij y^h'^j'A^ki) - 11 7T T-Trk+l,k+2{Ak, Aj) 11 TT — T-Trk,k+\{Ak, Aj) 

k=i tt2lAfc, Xj) a2{Xk, Xj) 

^n^#4T^^,fc+i(A.,A0 (88) 
k=i «2(Afc, Xi) 

Before discussing the results, we should note that the above expressions for multi-particle 
states indeed reproduce our previous findings for the two-particle (after considering appendix 
C) and the one-particle states. Now, from equations (82-84), it is direct to read of the n-particle 
eigenvalue expression, namely 

A(A.{A,}) = [.,(A)l-n'^#4 + MA)]-n-.^2|^ 

j=i Oig[Aj,A) -^-^ a7[A,Aj) 

-MXrfl-^^^^^^'\XAXl}) (89) 
j=i Q^glA, Aj) 

Following the same arguments given for the two-particle state, and in particular the dis- 
cussion presented at the end of appendix C, we easily derive that the unwanted terms vanish 
provided the rapidities satisfy the following Bethe Ansatz equations 

ui{Xi 



UJ2{Xi) 



A«(A = A,,{A,}), t = l,...,n (90) 
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Once again, the final results have been expressed in terms of the underlying auxiliary 
problem, which for a general multi-particle state is defined by 

where the inhomogeneous transfer matrix T^^\\, {Aj}) is 

T^'\K {X^})tl::Z = rl\7M, ^iK::iK A^) . . . f^riK An) (92) 

As we have commented before these results are direct extensions of those obtained for the 
two-particle state. We see that the Bethe Ansatz equations and the eigenvalues still depend 
on an additional auxiliary eigenvalue problem. In the language of condensed matter we would 
say that so far we managed to solve the "charge" degrees of freedom but still remains the 
diagonalization of the "spin" sector. As we shall see next the "spin" problem can also be 
solved in terms of the algebraic Bethe Ansatz approach. 



4.2 The eigenvalues and the nested Bethe Ansatz 

The task of this section is the diagonalization of the auxiliary transfer matrix T*-^''(A, {Aj}). 
For this purpose we have to set up another Bethe Ansatz which will result in "nested" Bethe 
Ansatz equations for the rapidities we began with. This problem, however, is equivalent to 
the solution of the 6-vertex model in presence of inhomogeneities and it has been extensively 
discussed in the literature (see e.g. refs. ^ 0)- Therefore we will only sketch the main 



steps of the solution for sake of completeness. First we write the transfer matrix T^^\X, {Xj}) 
as the trace of the following monodromy matrix 

TW(A, {A,}) = /:J\,„(A, Xn)C%,^_,{X, Xn-i) . . . C%,^{X, Ai) (93) 

where A^^^ is the two-dimensional "spin" auxiliary space. The Lax operator £^(\)^.(A, Aj) is 
related to the auxiliary matrix f(A, Xj) by a permutation on the x space and its matrix 
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elements are 

/I 0\ 

b{\,\j) 

a{\,Xj) b{X,\j) 

Vo 1/ 

We now go ahead applying the ABCD algebraic Bethe Ansatz framework 
inhomogeneous transfer matrix. Writing the monodromy matrix as 

A«(A,{A,}) fi«(A,{A,}) 
C«(A,{A,}) /^«(A,{A,}) 

and taking as the reference state the vector 





r«(A,{A,}) 



(94) 



for an 



(95) 



o™) = n' 



(96) 



we find the following relations 



AW(A,{A,})|o«) = |o«) 

n 

D«(A,{A,})|0«) = 
C«(A,{A,})|0W) = 



0(1) 



(97) 



The field B'^^\X, {^j}) plays the role of a creation operator over the reference state. To get 
its commutation rules we solve the Yang-Baxter algebra for the monodromy matrix T^^^ (A, {A^}) 
using as intertwiner the auxiliary matrix (26). This yields the following relations 

A«(A,{A,})i?«(/.,{A,}) = ^^5«(/.,{A,})A«(A,{A,}) 



D«(A,{A,})5«(/.,{A,}) 



i?«(/.,{A,}),i?«(A,{A,}) 



.g^i.«(A,{A.})^a)(MA.}) 
^ -i?«(/i,{A,})D«(A,{A,}) 



6(A,/i) 



(9^ 
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Next we have to make an Ansatz for the eigenstates of T^^\\, {Aj}). This is the "spin" 

m 

part of the multi-particle states and it is given by the product Y[B^^\fxi, lo^^^) whose 

1=1 

components are precisely identified with the coefficients jF'^" - "i. With the help of commutation 
rules (98) we are able to carry on the operators A^^\X, {Aj}) + D^^\X, {Aj}) through all the 
creation fields B^^\fj.i, {Aj}) leading us to the following result for the auxiliary eigenvalue 



A"'(A.{A,}.M) = n 



provided the numbers {/i;} satisfy the additional restriction 



1 = 1 m 



(99) 



(100) 



Finally, we use the auxiliary eigenvalue expression to rewrite our previous results for the 
eigenvalues and Bethe Ansatz equations of the "covering" vertex model. Substituting the 
expression (99) in equations (89,90) and using the second relation of equation (28) we obtain 
that the eigenvalue is 

- (X + [-3(A)]^ n -^^^TTTT 

a9(Aj,A) .^^ a7(A,Aj) 

1 " a.iXAi):^ 1 



A{X,{X,},M) 



[^i(A)]^n^"'^^'''^^ 



a9(A, Aj) f-J^ b{ni, A) 



"7(A, Aj) fj^ b{X,iJ,i) 



while the Bethe Ansatz equations for the rapidities {Aj} becomes 

L rn. -1 



UJ2{Xj 



n 



(101) 



(102) 



1=1 Aj) 

Now we are almost ready to make a comparison with the Lieb's and Wu's results 0. First 
we introduce a new set of variables z±{Xj) defined by 

6(A, 



z+{Xj 



(103) 



6(A,) a(A,) 
Considering this definition and taking into account the transformation (32) as well as the 
identities (29,30), we are able to rewrite the expression for the eigenvalue as 

■l + 2;_(A,)/^+(A)- 



-0"A(A,{^±(A,)},{^,}) 



l-^_(A,)/^-(A) 
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+[^3(A)]^ n 



6(A) 



l + z_(A,)z_(A) 



-MX)]' n 



A «(A) 
6(A) 



f=i 



l-z.{Xj)z4X) 

' l + z.{X,)/z+{X) 
1-^_(A,)A-(A), 



^ z4X)-l/z4X)-fii + U/2 
l\z4X)-l/z4X)-fli-U/2 



l\ «(A) 



l + ^-(A,)^-(A) 
l-z_(A,)^+(A) 



fT lA+(A)-.--+(A)-/^/-[//2 
y,l/z+(A)-z+(A)-/i, + C//2 

(104) 



and the nested Bethe Ansatz equations are now given by 



[^-(A,)]^ 



n 



z4Xj)-l/z4Xj)-ili + U/2 



t\z4Xj)-l/z4Xj)-fci-U/2' 



^ z4X,)-l/z^{X,)-f,i-U/2 ^ ^ I1,-I1, + U ^ l-l,...,m 



(105) 



J-J-^ z_ (Xj) - 1/z^ (Xj) - fli + U/2 ^J-^ jli-pLk-U' 

Prom the above expressions we note that function A{X, {z±{Xj)} , {fli}) is analytic in A. This 
happens because the condition of having zero residues on both direct Z-{Xj) and "crossed" 
z+{Xj) channels is clearly fulfilled by the nested Bethe Ansatz equations. The next step is to 
expand the logarithm of the eigenvalue A{X, {z±{Xj)} , {fli}) in powers of A and up to second 
order in the expansion we find 



ln[AiX,{z^{Xj)},{f,i})] 



-n 



+ ^ln[^_(A,)] + 



+A 



+A^ 



U 



^[^_(A,) + lA_(A,)] + -(L-2n) 



3=1 
n 



y: [z'_{Xj) - i/zi{x,) - c/[^_(A,) - i/^_(A,)]; 



3^ 



+C»(A^) 



-2L 
(106) 



The C(A) term parametrizes the spectrum of the Hubbard Hamiltonian and to recover 
the Lieb's and Wu's results we just have to reexpress the variable z_{Xj) in terms of the hole 
momenta kj by 



z_(A,) 



(107) 
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Considering this relation, the eigenenergies of the Hubbard model are 

En{L) = ^^^"^""^ + j2 2 cos(A;,) (108) 
and the momenta kj satisfy the following Bethe Ansatz equations 



e 



iLk 



n 



sin(fcj) — fii — iU / 4 



\ sm{kj) - fii + iU/A' 



y sin(A;,) -Jii- lU/A Jii-JJik + iU/2 ' 



m (109) 



where we also used jli = 2ifii to bring our equations in the Lieb's and Wu's form. A careful 
reader might note that the above Bethe Ansatz equations have an extra minus factor in front of 
the coupling U in comparison to the original ones. This is because we are using the language of 
holes instead of particles and this means that the integers n and m are the total number of holes 
and the number of holes with spin up, respectively. It is well known that via a particle-hole 
transformation the kinetic term of the Hamiltonian gets an extra minus sign, which changes 
the sign of factor U/t entering in the Bethe Ansatz equations. Similar reasoning can be carried 
out for others conserved charges. For example, the first non-trivial current commuting with 

is 



the Hamiltonian IIT2|, 113 



i=i 

+U{c]^Cj+^ - c]+i^9|)(nj-+u + Tiji - 1) + [t^i] (110) 
and from equation (106) it follows that the spectrum (modulo a constant) of this charge is 

n 

Ei{L) = 2iJ2[sH^kj) - Usm{kj)] (111) 
i=i 

We would like to close this section commenting on the construction of the eigenvectors in 
the terms of the "dual" field B*{\). The equivalence between the commutation rules for the 
fields B{X) and B*{X) allow us to follow straightforwardly the whole construction of section 4.1 
and it is not difficult to derive formula for the "dual" eigenvectors $*„(Ai, . . . , A„). Formally, 
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we can apply the "dual" transformation in expression (81). This leads us to following "dual" 
recurrence relation 

#„(Ai,...,A„) = #(Ai)®$Vi(A2,...,A„)+X:^"*°^^''^'^ 



i=2 



X ® F(Ai)$'*„_2(A2, . . . , A,_i, A,+i, . . . , A„)Z^(A,)] (112) 

We expect that the corresponding eigenvalues A* {X, {Xj} , {fli}) should also be related to 
A(A, {Aj}, {fti}) in some way. This is indeed the case if we shift all the rapidities around the 
"crossing" point tt/2, and the relation we found is 

A*(7r/2 - A, {n/2 - A,}, {n/2 - fn}) = (-1)"A(A, {A,}, {fn}) (113) 

With this we complete our analysis of the graded eigenvalue problem and in the next section 
we shall discuss some other complementary results which can be obtained within the ABCDF 
formalism. 



5 Complementary results 

In this section we shall first consider the solution of the coupled spin model with twisted 
boundary conditions. This allow us to illustrate the difference between the Hubbard and the 
coupled spin models from the viewpoint of their Bethe Ansatz solution. Next we consider 
the well known SU{2) symmetries of the Hubbard model P^, BBl 15^. We will show that 



the eigenvectors (81) are highest weights of both the SU{2) Lie algebra of rotations and the 
ry-paring SU{2) symmetry. Thus we are able to recover the results by Essler, Korepin and 
Schoutens [B^ from an algebraic point of view. 
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5.1 Twisted boundary conditions 

We begin recalling that twisted boundary conditions are in general associated to certain gauge 
invariances of the Yang-Baxter algebra. The integrability condition (1) is still valid when 
'^AiW G_A^C_AiW provided the gauge matrix G_4 satisfies p5| 



[R{X,fi),G^®G^]=0 (114) 

This means that a vertex model defined by the transfer matrix Tg(A) = TrXTci^) whose 
monodromy matrix is 

Tg{\) = Gj,Cj,l{\)Cm-i{\) ■ ■ ■ Cm{\) (115) 

still remains integrable. One way of seeing the connection to twisted boundary conditions 
is, for example, to derive the quantum Hamiltonian Hq commuting with the transfer matrix 
Tg'(A). To this end we assume that the Lax operator is regular at some value of the spectral 
parameter, say C^AiiS^) = Pai where P^Ai is the ® permutation operator. Then the local 
quantum Hamiltonian Hq is [|], 

HG = T^'m'^{^) (116) 
where symbol ' stands for the derivative on A. By using the permutation properties 

dMPA,=PAjd,,, PMdAj=d-,,PM, Trj,[GAPAa]=G^ (117) 



and after few algebraic manipulations we derive that (see e.g. ref. |^ 



L-l 



H=Y. + GJ^Hl^^Gl (118) 

i=l 

where hij = [PC'{0)]ij and we assumed that is invertible. The last term in the Hamiltonian 
(118) refiects the presence of non-trivial boundary conditions. In the context of the coupled 
spin model (7), it is straightforward to see that we get twisted boundary conditions by taking 
the following gauge 

^•^ = I )^[ e^<^^/2 ) ^^^^^ 
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Clearly, such gauge matrix fulfill the integrability condition (114). In order to diagonahze 
Ta{X) we only need to introduce few modifications on the formalism developed in the previous 
sections. It is fundamental that this gauge does not spoil the triangular form of the monodromy 
7g(A) when it acts on the ferromagnetic reference |0). The diagonal operators of ^(A), 
however, pick up extra phase factors and now we have the following relations 

B{X) |0) = e-^('^i+'^^)/2[^^(^)]L ^^^^ _ e^('^i+'^2)/2[^^(^)]L 
An (A) |0) = e-'^'^'-'^'^/^MX)]'' |0) A22(A) |0) = e^^^^-^^^/^lc^alA)]^ |0) 

(120) 

The next step is to solve the commutation rules in the standard Yang-Baxter formalism, 
since we are considering the coupled spin model. These commutation rules have basically the 
same structure of those worked out for the graded case, apart from few signs and imaginary 
factors. We have collected them in appendix B and we note that the corresponding 6-vertex 
auxiliary matrix has now an extra sign in the amplitude b{X,fi). Therefore, the nested part 
always gets twisted, emphasizing the difference between the Hubbard and the coupled spin 
models for closed boundary conditions. Since now the basic ingredients have been set up we 
can follow closely the steps of sections 3 and 4. Here we are interested in the eigenvalues of 
the twisted model and now we begin to summarize our final findings. Taking into account the 
relations (120) and the commutation rules B.12-B.22 we derive that the eigenvalues of transfer 
matrix To (A) is 

Ae(A,{A,}) = e-('^^+'^^)/2[a;i(A)]^n-^4T^ + ^'^^'^^'^^'[^3(A)]^n ""'^^'^'^ 



j=i oig(Xj,X) .^^ a7(X,Xj) 



_^ a^^g.j^ (121) 

j=l Oig{X,Xj) 



where the variables {Xj} satisfy the Bethe Ansatz equations 



ui{Xj 



L 



, -(-l)V(^^+^^)/2A«(A = A,,{A,}), i = l,...,n (122) 

It turns out that the auxiliary problem gets also an extra modification besides the sign on 
amplitude 6(A, //). The auxiliary problem absorbs the twisting on the diagonal fields Aii(A) 
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and ^22 (A) and now function Aq\X, {Xj}) is the eigenvalue of the following auxiliary transfer 
matrix 



(123) 



where the Lax operator >C^j(A, Xj) and the matrix G]^' are given by 
/I 0\ 



' ^-i{4>i-4>2)/2 g 

gi(01-</.2)/2 



(124) 



-b{X,Xj) a{X,Xj) 
a(A,Aj) -6(A,Aj) 

yo i; 

The solution of this auxiliary problem is once again standard. Following the lines of section 
4.2 we find that the auxiliary eigenvalue expression is 

A« (A, {A,}, {/.,}) = e-(^-^^)/^ n ^T^T + e^^^-^^)/^ R "KA, A,) U -^-^ (125) 

Collecting these results altogether and substituting the variables z±{Xj) and jli we find that 
the eigenvalues of 7g(A) can be written as 

6(A)1 + ._(A,)/,:^(A) 



Ag(A,{;.±(A,)},{AJ) 



(^^^^^)/^[..3(A)]^n^i 



+e' 



j^XA)l-z_(A,)A_(A) 

+ --(A,)--(A) 
-z4Xj)z+{X) 



X 



n 



\z4X)-l/z4X)-jli-U/2 

6(A) 1 + z_(A,)^-(A) ™ l/^+(A) - ^+(A) - - U/2 



n 



a(A) 1 - ^-(A,)^+(A) t=\ l/^+(A) - ^+(A) - /i, + t//2 



(126) 



while the nested Bethe Ansatz equation are given by 



-l)'"+"e-'^2[;2_(Afe)]^ 



z_{Xk)-l/z4Xk)-fli-U/2 
j^J--, z4Xk) - l/z4Xk) + U/2 



^ z4X,)-l/z4X,)-flj + U/2 ^ A;=l,...,n 



n 



z_(Afe) - l/z4Xk) - flj - U/2 

.,. , . in u, - Tu + 

:, I = 1, . . . ,m 



=1 fii-H-U 



(127) 
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In order to get the results for the Hubbard model with twisted boundary conditions we 
substitute the angles (10) in the above expressions. We should also remember that we are 
using the language of holes and therefore the integers n and m are identified with the total 
number of holes A^'* and the number of holes with spin up N^, respectively. This cancels extra 
phase factors in the Bethe Ansatz equations (127) and we recover the known set of nonlinear 
equations parametrizing the spectrum of the twisted Hubbard model B^, HO]. Let us close 



this discussion by mentioning a possible application of these twisted Bethe Ansatz results. 
Consider the Hubbard model perturbed by a particle current term (see e.g ref. |^) with 
periodic boundary conditions. This model is described by the Hamiltonian 

L 

H,iU, K) = HiU, 0t = 0, 0i = 0) - E - cLq+1^) (128) 

i=l o-=± 

In the spin language, this perturbation is a Dzyaloshinsky-Moriya interaction in the az- 
imuthal direction, playing the role of a "vertical" magnetic field. Similar to what happens 
in the spin case |^2|, the fermionic current perturbation can be gauged away by using the 



canonical transformation [ED 



. (2k-3)(l> 

Cka 2 Cfca , tan(0) = Ac (129) 



allowing us to derive the relation 



Hc{U, A,) = Vl + Agg( , ^ , 0T = 0^, 0i = <i>L) (130) 

1 + A^ 



Thus, the spectrum of Hc{U, Ac) is related to that of the Hubbard model with certain twisted 
boundary conditions and renormalized coupling. Similar reasoning also works if we add a spin 

L 

current term Jg = —iXg E E — cfo-Cj+io-) to the Hubbard Hamiltonian. In this 

i=l a=± 

case, after performing the transformation 

. (2fc-3)0 . (2fc-3).j> 

Cfct e* 2 , Cki ^ e ' 2 (131) 

we find the that the Hamiltonian Hs{U, A^) of the Hubbard model perturbed by the spin current 
satisfies 

U 



Hs{U, Xs) = + mi I ^ , 0T = <PL, 0^ = -0L) (132) 

'1 + \i 
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Before closing this section we would like to comment on possible extensions of the results we 
have obtained so far. First it is possible to diagonalize a two-parameter family of vertex models 
whose Lax operator is C^^'^\X) = PR{X, 9q) O, E^. Its Bethe Ansatz solution follows directly 



from the results of this section, since the main change is only concerned with the action of the 
fields on the reference state. It turns out that now the bare pseudomomenta (left-hand side of 
first equation (127)) depends on the variable 6*0 as [i^f^^^^]'^- Also, the whole formalism can 
be extended to treat the Hubbard model in the presence of chemical potential |Q. Finally, 
for further results on twisted boundary conditions see for instance ref. p4| . 

5.2 SU{2) symmetries 

In this subsection we investigate the highest weights properties of the eigenvectors constructed 
in section 4, with respect to the two SU{2) symmetries of the Hubbard model Few years 



ago, Essler Korepin and Schoutens |^ have shown that certain "regular" states obtained from 
the coordinate Bethe Ansatz wave function are highest weight states of both the SU{2) algebra 
of rotations and r^-pairing SU (2) symmetry. The idea here is to explore the algebraic machinery 
we developed in the previous section to study this problem from an algebraic perspective, in 



close analogy with the discussion by Takhtajan and Faddeev|[45|| for the Heisenberg model. For 
this purpose we will use the results of Gohmann and Murakami who recently showed that 
the graded monodromy matrix indeed commutes with these two SU(2) Lie algebras. More 
precisely, following the notation of ref. we have 



[r(A), = -[r(A), « = +,-, ^ (133) 

and 

[r{X),r]Q = -[T{X),J2U, a = +,-,z (134) 

where the subscripts Q and A emphasize in which space, quantum or auxiliary, the commuta- 
tors are taken, respectively. The SU(2) generators of rotations 5" and those of the //-pairing 
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symmetry 77" are defined by [p7 



and 



5+ = -EsV.U , ^" = -EsV.T > ^^ = E(^.T-%i) (135) 
j=l j=l j=l 



^"^ = E(-1)''^Mt4 ' r = El-iy-^'c.iC.T , r/^ = EKT+^.i-l) (136) 
j=i j,=i j,=i 



while the matrices X)" and J2 a,re ||3^ 



^+ = a+®a- , J2 =(y~®(r+ , = ^((T^®/-/®a^) (137) 

Y^^ = a+®a+ , ^~=a-®a- , = + J®a") (138) 

Let us begin by considering the r^-pairing symmetry. The identity (134) enables us to 
compute the commutators of the creation fields B{X) and F{X) with the SU{2) //-pairing 
generators. For the component 77^ we find 

[r/^5(A)] = -5(A) , [r/^F(A)] = -2F(A) (139) 

while for rj'^ we have 

[V+,B{\)] = ~C*{\) , [v\F{\)]=B{\)~D{\) (140) 

We see that formula (139) corroborates the physical interpretation we have proposed for 
the creation fields -B(A) and F{\), i.e. that they create a single and a doubly occupied hole on 
the full band pseudovacuum. For example, from this equation it is straightforward to derive 

ri' |$„(Ai, ■ ■ ■ , A„)) = {L-n) |$„(Ai, ■ ■ ■ , A„)) (141) 

where we used the property r/^ |0) = L |0). 

We note that the above result is valid for arbitrary values of the rapidities. However, this 
is no longer true when we consider the annihilation property of the raising operator r/"*". In 
what follows we shall show that 

r/+|$„(Ai,---,A„)) = (142) 
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provided the rapidities {A^} satisfy the Bethe Ansatz equations derived in section 4. 

To verify the above annihilation property it is instructive first to study the case of few 
particles over the reference state and afterwards use mathematical induction for the general 

case. From equation (136) this is clearly correct for the reference state. For the one-particle 
state, by using the first commutator (140), it is easy to show that 

ri+ |$i(Ai)) = v^B{\,)J |0) = [7?+, B{\,)]J |0) = |0) = (143) 

The Bethe Ansatz restrictions start to emerge in the two-particle state analysis. For this 
state the commutators (140) produce 

r;+|$2(Ai,A2)) = B{\i)(E)r]+B{X2).f\0)-C*{Xi)^B{X2).f\0) 

+ '-^i^^mXi)-D{X,)]B(X2)lf\0) (144) 

tt7(Al, A2) 

The first term in the above equation vanishes by the same arguments used in the one- 
particle state analysis. To simplify the second term we use commutation rule (B.3) and finally 
the third term is easily estimated from the diagonal relation (23). Putting these simplifications 
together we find 

r7+|$2(Ai,A2)) = ^^i^^^^fK(A>i(A2)]^-[w2(Ai)w2(A2)]^le-^|0) 



a7(Ai, A2) 
iaio(Ai. A2) 



[wi(Ai)wi(A2)] 



L 



avi^i, A2) 

-k2(Ai)^2(A2)]^A('nA = Ai, {A^})A«(A = A2, {A,})] If |0) 

= (145) 

where in the second fine we used the following two-particle identity, A(^)(A = Ai, {A;})A(^)(A = 
A2, {A;}) = 1. Clearly, the term in brackets vanishes due to the Bethe Ansatz equations (90). 

Next we consider the three-particle state. We shall see that a general pattern in the analysis 
begins to emerge here. After using the commutator relations (140) we have 

77+|$3(Ai,A2,A3)) = S(Ai)®r7+$2(A2,A3).J?|0)-C*(Ai)®$2(A2,A3).J^|0) 

+i[B{Xi) - D(Ai)] ® B(A3)5(A2)(?r (Ai, A2, A3)./ 10) 
+i[B{X^)-D{X^)]^B{X2)B{Xs)gi'\X^,X2,X^)■f\0) (146) 
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The first term is computable directly from the first line of equation (145), after making the 
replacements Ai — > A2 and A2 — > A3. The third and fourth terms are estimated with the help 
of commutations rules (35-36). The simplifications for the second term is more complicated 
since it involves the knowledge of an extra commutation rule, besides relation (B.3), between 
the fields C*{X) and -F(//). This relation is given by 

^U{X)D{t,) + ^^^DiX) (147) 



q;7(A,/x) q;7(A,/x)' 

Collecting all the pieces together is remarkable to see that many terms have opposite signs 
and thus they are trivially canceled out. However, there is a non-trivial simplification yet to 
be carried out. This is related to the terms proportional to [^.{B*{Xi) (8) /)] and they vanish 
thanks to the following identity 

aio{x,z) ag{x,y) ^ a^jy, z) aio{x,y) ^ aio{x, z) a2{z,y) ^ ^ 
arix,z) arix,y) ag{y,z) a^ix^y) a-jix^z) ag{z,y) 

After these simplifications, the remaining terms are only proportional to <S> B{Xj)] and 
they can be compactly written in the following way 

77+ |$3(Ai, A2, A3)) = E ■ ■ ■ , Ai, . . . , A,-, . . . , Xs)]Q\f{Xi, Xj; {A,}).J? |0) (149) 

i=2 1=1 

The first term (5i2^(Ai, A2; {Ajt}) is easily figured out because it has only two main con- 
tributions coming from the second and the third terms of equation (146). The other two 

are obtained from this term via consecutive permutation of rapidities through the exchange 
property (78). The expressions for these coefficients are 

QlfiXi, A,; {Xk}) = [K(AOti;i(A,)]^ - [w2(A0w2(A,)]^A«(A = A^, {Afc})A«(A = A,-, {A^})^ 

a^X^ Jl ^4^^#44o;f(Ai,A,;{A,}) (150) 
ariXiAj) l\ iag{Xi,Xk)iag{Xj,Xk) ^ ''^ ^ ' 

and they vanish again as a consequence of the Bethe Ansatz equations (90). 



42 



Now using mathematical induction it is possible to write the action of the raising operator 
on a general n-particle state as 

ri+ |$„(Ai, . . . , A„)) = ^ Y^ii® 4-2(Ai, . . . , A^, . . . , A,-, . . . , A„)]gJ;)(A,, A,-; {Xk})J |0) (151) 

i=2 1=1 

As before, it is convenient first to compute the simplest coefficient Q12 {^i^ ^2] {^k}) and 
then take advantage of the permutation property (78) to obtain the remaining ones. For this 
term we have just two contributions coming from 

I 5(Ai)4-2(A3, . . . , Xn)B{X2)g[''\Xi, ...,Xn)-^ |0) (152) 

and 

II := -(7*(Ai) ® 4-i(A2, . . . , A„)].J? |0) (153) 

We compute the first part by carrying the scalar operator -B(Ai) through the vector 
$„_2(A3, . . . , A„) keeping only the "wanted terms" proportional to B{Xi). This is very similar 
to what we did in appendix D and we find 

T I /\ \ /\ mL^"1o(^1, ^2) -A- ?'tt2(Afc,Ai)?:a2(Afc,A2).-* - , M'^ln\nc^^ 

7 := [w;i(Ai)m;i(A2)] [[ — — — — -[^ ® $„-2(A3, • • • , A„)].J^ |0) (154) 

Qi7(Ai,A2) Q;9(Ajfc,Ai) Q;9(Ajfc,A2j 

fc^l,2 

The second part is more involving since we have to carry two operators of type ^(A) through 
vector $n_2(A3, . . . , A„). This means that we have to compute the expression 

J J - '"f,^'\^'^ e./3ia5.(Ai)i/36.(A2)[^n-2(A3, ■ ■ ■ , A„)],3...,„^''"-^^ |0) (155) 

which after some algebra can be compacted back as 

iQ;io(Ai, A2) n s s,^ -A- ai(Ai,Afc) ai(A2,Afc) 
^^■^ (\ \ \ [^2(Ai)w2(Ai)] II . . ^ . . ^ 

Q;7(Ai,A2j ZQ!9(Ai, Afe) ZQ!9(A2, Afej 

fe^l,2 

x[f ® C-2(A3, . . . , A„)]«,...„„[T(^)(A = Ai, {Xi})T^'\X = A2, {Az})]^V±„-^'"-'^ |0) (156) 

Finally, putting together expressions (154) and (156) and also using the auxiliary eigenvalue 
definition (91) we find 

gS?(Ai, A2; {Xk}) = [[w^{X,)w^{X,)]'^ - [w2iX^)w2{X2)]'' A^'\X = X,, {Xk})A^'\X = A2, {Xk}) 

iaio(Ai, A2) -pj- Q:i(Ai,Afc) Q;i(A2,Afc) (157) 
tt7(Ai,A2) l}^ iag{Xi, Xk) iag{X2, Xk) 
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which once again vanishes due to the Bethe Ansatz equations. All the other coefficients are 
obtained by permuting the rapidities and by taking into account the exchange property (78), 

^ (2) 

and as a result they get an extra multiplicative "ordering" factor O^y (A;, A^; {Xk})- Since the 
Bethe Ansatz equations are invariant under indices relabeling, they vanish too. This completes 
the proof that the eigenvectors (81) are highest weight states of the r^-pairing symmetry. 

Next we turn to examine the highest weight property of the SU (2) algebra of rotations. Now 
the commutators of the creation fields with the SU{2) generators are obtained from equation 
(133). For the component we find 

[^^i?l(A)] = 5l(A) , [S\B2{X)] = -B,{X) , [S\F{X)] = (158) 

and for we have 

[5+,i?i(A)] = 0, [S^,B2{X)]=B,{X), [^+,F(A)] = (159) 

First of all, it is not difficult to see that eigenvector (81) will be hardly annihilated by the 
raising operator unless further restriction are assumed. To illustrate this fact in a simple 
example let us consider the one-particle state. By using the commutators (159) we find 

5+|$i(Ai)) = 5i(Ai)^2|0) (160) 

where we used that |0) = 0. Therefore, to assure the highest weight property for the one- 
particle state we must set JF^ = 0. This is an example of what was called "regular" Bethe states 
in ref. |3^, and in general these states are obtained by projecting out the negative sectors of 



the magnetization operator S^. This later condition is easily implemented for the eigenvector 
(81) if one uses the commutators (158). 

To see how this works in practice let us consider the two-particle state. In this case it is 
obvious that we have to set JF^^ = , and after that we find 



|$2(Ai, As)),,,,;,, = S+ [5i(Ai)5i(A2)^^^ + fii(Ai)52(A2)^'^ + B2{X,)B^{X2)J'''+ 

|0) 

E^^']5i(Ai)5i(A2)|0) (161) 



+,^l!j^F{X,)B{X,)lf 
a7(Ai,A2) 
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Table 1: The "regular" multi-particle states properties up to n = 4. 
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^1122 _ ^ J7II22 _ ^ J7I212 _ ^ J72112 _ Q 
P P P P 



where the sum is over permutations on the indices of the coefficient ^"^"i. In this case it 
is straightforward to verify that this sum indeed vanishes by directly solving the auxiliary 
eigenvalue problem (56). The deeper reason behind this fact, however, is that the vanishing 
of such sum is precisely related to the highest weight property of the Bethe wave functions 
of the XXX Heisenberg model with two sites. We should recall here that the components of 
this wave function are identified with the coefficients ^"201 Yxoth this discussion, it becomes 
evident that the whole procedure can be applied to any multi-particle state. As an example, 
in table 1 we summarize our findings up to the four-particle state 

The columns of table 1 refer to the particle number, magnetization values and the sufficient 
vanishing condition for S'^ annihilate the "regular" part of eigenvector (81), respectively. In 
the sum the symbol a means that the a-th element is maintain fixed under permutations. 
The generalization to multi-particle state is done by induction and the sufficient vanishing 
conditions are made of the many possible permutation over the coefficients J^°-^-°-^ having 
posit ve magnetization. As before, these conditions are fulfilled as a consequence of the highest 
weight property of the Bethe states of the XXX Heisenberg spin chain in a lattice with size n. 
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Since this later point has been well explained by Essler, Korepin and Schoutens there is no 
need to proceed with details, and thus we conclude our proof that 5*+ |<l'„(Ai, ■ ■ ■ , ^n))reguiar ~ ^ 
here. 

Finally, we remark that similar properties can be also verified for the "dual" eigenvec- 
tor. The only difference is that now the "regular" states are defined by projecting out the 
positive sector of the magnetization. At this level, the eigenvector and its "dual" becomes 
complementary eigenstates. 



6 The ABCDF framework for the Bariev model 



The purpose of this section is to illustrate that the ABCDF framework developed in the 
previous sections is by no means only applicable to the Hubbard model. In order to show 
that, we consider a second interesting model of interacting XY chains whose corresponding 
i?-matrix also does not have the difference property. The model was originally formulated by 
Bariev |^ and its one-dimensional Hamiltonian is 



(162) 



i=l 



where is a coupling constant. In the language of fermions V plays the role of a bond-charge 
interaction and Hamiltonian (162) resembles the model of hole superconductivity proposed by 



Hirsch |^ 



In the context of the quantum inverse scattering method this model has recently been 



investigated by Zhou |^ and Shiroishi and Wadati |^0[ who found two distincts covering 
vertex models for the Bariev Hamiltonian. In this section we apply the ABCDF formalism 
for the former solution^. In this case, the proposed Lax operator was |^| 



"^Part of our results were first announced in ref. Esl. See also ref. IBll 



'(1), 



'(2), 



(163) 
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where 



^S(A) = ^(l + a|a^) + ^(l-a|a^)exp(/?r+r^) + (a+a^ + a7a+)Vl + A2exp(2/3r+r^) (164) 



A 



and 



>C3](A) = ^(l + r;r^) + ^(l-r;r^)exp(/?a>^) + (r;r^+rrr+)Vl + A2exp(2/?a>^) (165) 
The relation between the parameter jS and the couphng constant V is determinated by 

cosh(/3/2) 



A 



computing the expression P-^C^^\\) on A = 0. After performing the rescahng A 



we found 



h 



V 



1-V ^'''^ 
The i?-matrix solving the Yang-Baxter algebra for this choice of Lax operator was also 
found by Zhou. Its explicit 16 x 16 form is E7 
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(167) 



where the fifteen non-null Boltzmann weights pj{X,fi), j = 1, . . . , 15 have been collected in 
appendix E. We remark that we have verified that this i2-matrix indeed satisfies the Yang- 
Baxter equation (6). 

We note that the structure of such i?-matrix is very similar to that found for the Hubbard 
model and consequently one could easily guess that the ABCDF formalism should work for 
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this embedding as well. It is not difficult to adapt the main steps of section 3 in order to 
obtain the appropriate commutation rules for such classical vertex analog of the Bariev model. 
The most important commutation rules have been summarized in appendix E. The interesting 
feature here is the structure which comes up for both the "exclusion" vector and the auxiliary 
r-matrix. We found that they are given by 

/I 0\ 

a(^)(A,^) 



(0 1 l/h 0), f(^)(A,/^) 



a[''\X,i2) 

Vo ly 



(168) 



where the weights a^^\\, fi) and b^^\\, /x) are 
A(l - 



(169) 



A — /i^/i A — K^ii A — K^n 

From equation (169), it is easily recognizable that the auxiliary r-matrix has the structure 
of an asymmetrical and anisotropic 6-vertex model because the parametrization leading to 
the difference property for f^^\\, fx) is now standard, namely A = exp(iA;). In this case the 
hidden symmetry is of Hecke type because such auxiliary r-matrix can be produced as a result 
of Baxterization of the Hecke algebra (see e.g. ref. |^2|). We recall here that this later 
symmetry was first noted by Hikami and Murakami by exploiting the continuum limit of the 



Bariev Hamiltonian ||^. Interesting enough, we note that the "exclusion" statistics for "spins" 
degrees of freedom seems to be of anyonic type with a phase jS which depends on the strength 
of the coupling constant V (see equation (166)). It remains to be seen if this feature will 
also be manifested in physical quantities computable by Bethe Ansatz methods such as in the 
low temperature behaviour of the free energy (conformal limit) and in the scattering of the 
elementary excitations. 

Let us now discuss the construction of the eigenvalues and the eigenvectors for this classical 
analog of Bariev model. It turns out that such formulation goes fairly parallel to the one already 
presented in section 4 and in appendix D. For this reason we shall avoid unnecessary repetition, 
and from now on we concentrate our attention only to the basic points. We start directly with 
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the two-particle state analysis since it has already proved to contain sufficient information about 
the main steps entering in the relevant computations. Afterwards, generalization to multi- 
particle states is made following similar discussion presented in appendix D. Our previous 

experience with the Hubbard model suggests us to begin with a symmetrized two-particle 
vector. As before, the main trick is to look at the commutation rule between the two creation 

— * 

fields of type B{X). Prom equation (E.17) it is not difficult to guess that such vector is 



We go ahead computing the action of the diagonal fields on the two-particle state Ansatz. 
Here we shall use fully the permutation property of the eigenvector, specially the simplifications 

mentioned at the end of appendix C. Considering the commutations rules of appendix E and 
following the calculations of section 4, we find that the expressions for the action of the diagonal 
fields on the two-particle state are 




(170) 



which is indeed the case thanks to the following identity 




(171) 




2 



i:|f^«(A,A,;{A,})) 



i=i 




(172) 




2 



- E[^^#A« (A = A„ {A,}) A,; {A,})) 



Ha, Ai, A2)[^'AiA2]^ Xj, Xi; {Xk})) (173) 
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-E|f^S'HA,A,;{A.})) 

+//f )(A, Ai, A2)[;iA2]^Ag)(A = A2), {A,}) \d^\x, A,-, A^; {A,})) 
+i/f )(A, Ai, A2)[/iAi]^A« (A = Ai),{A,}) (A, A,-, A,; {A,})) 

(174) 

where we used the relations 5(A) |0) = |0), A^aiX) = [hX]^ |0) and D{X) |0) = [A^ |0) which 
are determined by acting the Lax operator on the ferromagnetic pseudovacuum. As before, 
A^^^A, {A;}) is the eigenvalue of the auxiliary problem (56) whose r-matrix is now r^^\X,iJ,). 
Furthermore, the expressions for the unwanted terms are 



Q?{X,Xf,{Xi}) 



P2(Aj, A) -j^ pi(Afc, Aj) 
P3(Aj-,A) P3{Xk,Xj) 



S(A)®S(A,)]nCii(Afe,A,).:?|0) 



(175) 



1]P(A,A,;{AJ) 



(176) 



oJ'Ha, A,, A,;{A,})) = F(A)e^^)./|0) 

Finally, the functions H( (a;, 2;), / = 1, . . . , 4 are given by 

r(s)/^ _ P5(l/,^)P4(l/,a;) , pi(?/,x)p2(^,a:)pi2(|/,a;) 



^ ' ' P9{x,y)p9{y,z) pg{x,z)pg{x,y) 

Pi{x, y)p2{x, z)p5{x, y) p^{x, y)p2{x, y)p2{y, z 



pg{y, z)pg{y, x) psiy, x)p3{z, x)pg{y, x) 

P4(.r. U)P^>{U- z) p^{x. z)p2{,r. tj) 



Hi''\x,y,z) = [ 



(177) 



(178) 
(179) 



Paix, y)pz{x, z)p9{x, y) pg{x, y)pz{x, y)pz{y, z) 



][a''\y.z) + ''^-^] (180) 



h 



Hf\x,y,z) 



Pi{.r. -)p2(.'-- y)p-,{'''- z) _ p5{-c- '-)p2(.'-. --)P2(--- y) 

Psix, Z)p3{x, y)pg{x, Z) pg{x, z)p3{x, z)p3{z, y) 
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(181) 



In order to cancel out the unwanted terms it is sufficient to impose the following Bethe 
Ansatz restriction to the rapidities 

[A,/i]-^ = -Aij')(A = A„{A,}), 2 = 1,2. (182) 

since this condition eliminates automatically the first two kind of unwanted terms. Moreover, 
this helps us to gather the four unwanted terms proportional to F{X)^^^\^ which are finally 
vanished due to the identity 

i/J^) {x,y,z) + iff) {x,y,z) = fff ^ {x,y,z) + i/f ^ {x,y,z) (183) 

To obtain the two-particle eigenvalue we collect the wanted terms and by using the expres- 
sion for the Boltzmann weights (see appendix E) we find 

A™(A. {A.}) = n 'l^y + n + n ''''+^/'\ W(A. {A,}) (184) 

j=l - ^ j=l ^- j=l ^ ~ 

The generalization of these results for multi-particle states goes much along the lines dis- 
cussed in appendix D. We start constructing a symmetrized n-particle vector state which 
satisfies 

$(^)(Ai, . . . , A„ A,+i, ...,Xn) = ^^^HAi, . . . , A,4-i, A,-, . . . , AO.f(^)(A„ A,+i) (185) 
and after solving these constraints we have 

$<f\K . . . . A„) = B(A,) e\(A.. ...,K)-t ^4H4 n ^#44 

j=2 P9[M, ^j) fe^2 P9[^k, ^j) 

X 'i® F(Ai)CkA2, . . . , A,_i, A,+i, . . . , A„)i?(A,)] CVi(Afc. A,) (186) 

fc=2 

From the two-particle analysis it is not difficult to see what should be the expressions for the 
multi-particle eigenvalues and Bethe Ansatz equations. For example, the auxiliary eigenvalue 
expression is the same as given in equation (99), replacing 6(A,/i) by b^^\X,fj,). To make a 
comparison with the previous Bethe Ansatz results derived by Bariev EM it is convenient to 
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redefine the spectral parameter A, the rapidities {Aj} and the nesting variables {^j} |4^- Here 
we set 

A = hXj = e'''\ fij = e'^^ (187) 
In terms of these new rapidities, our final results for the eigenvalues are 

A{k, |A:,|, |A,|) - n ^ ^.^^^^^2 _ + ^/3/2) + ^'^^^^ z sin(fc/2 - h/2 + 

+ exp H(k - ^^)L] lfTl^2<tIl±Nl^im n sin(A,/2-A;/2 + z/3) 
+ exp Kfc z/^)LJ I n ^.^^^^^2 - A;/2 + ^f3/2) l\ sin(A,/2 - k/2) + 

" zcos(fc/2 + fc,/2-z/3/2) - sin(fc/2 - A,/2 + zP) \ 
l\ sm{k/2-h/2 + tP/2) sin(A;/2-A,/2) ' 

while the nested Bethe Ansatz equations for the rapidities {fcj} and {Aj} are 

exp(*,L) ^ -(-l)"n -'t/2"A/r'!!/2! - -l-" 

^^^^^ sm(fci/2 - Aj/2 - z/?/2) 

- sin(A,/2 - fc./2 - ^/3/2) ^ - sin(A,/2 - A,/2 - 
^ ^ iJfsin(A,/2-A;,/2 + ^/3/2) sin(A,/2 - Afc/2 + ' ^ ^ 

Finally, to obtain the eigenspectrum of the Hamiltonian (162) we expand the transfer matrix 
eigenvalues in power of the spectral parameter. Up to second order we have 

n -I n -1 

/n[A(A, {A,}, {/i,}] = E )^ + hX^Ch>^^ + jj-) 

+^i:[ijJ^r-ihX.?] + OiX') (190) 



Considering the 0{X) term of the above equation and remembering to perform the rescaling 

cosh{/3/2) 



A — » '^1(^/21 conclude that the eigenenergies of the Hamiltonian (162) are 



En = 2{1 + V) Y^cosih) (191) 

i=l 

We conclude remarking that this model can also be solved with twisted boundary conditions 
following precisely the same steps presented in section 5.1 . 
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7 Conclusions 



The main purpose of this paper was to apply the quantum inverse scattering program for the 
one-dimensional Hubbard model. We succeeded in developing a framework which allowed us 
to present an algebraic formulation for the Bethe states of the transfer matrix of the classical 



"covering" Hubbard model proposed earlier by Shastry [p^ |T3]. A hidden 6- vertex symmetry 
has been revealed, and it played a fundamental role in the solution of the transfer matrix 
eigenvalue problem. We have found the eigenvalues of the transfer matrix and showed that its 
eigenstates are highest weights states of both the rotational and the //-paring SU{2) symmetries. 



This later result corroborates the original proof given by Essler, Korepin and Schoutens |36 
in terms of coordinate wave functions. We have also discussed the algebraic solution of models 
with twisted boundary conditions and applied the results to the Hubbard model perturbed by 
charge and spin currents. 

The framework developed in this paper, the ABCDF formalism, is indeed suitable to solve 
a broad class of integrable systems. As an example, we solved, in section 6, the classical analog 
of the Bariev model by this method. There are also other models that fit in the ABCDF 
framework, such as the trigonometric vertex models based on the -B„, C„ , Dn , and A\^_^ 
algebras as well as certain related supersymmetric models . Interesting enough, the former 
models almost exhaust the Jimbo's and Bazhanov's list of Uq{G) R-matrices |5^, and only 
the -D^+i model appears to be not solvable within our framework. Anyhow, these examples 
suggest us that the ABCDF formalism is capable of solving integrable models having one less 
trivial conserved quantum number when compared to the An multi-state 6- vertex models with 
an equivalent Hilbert space. 

Finally, the possibility of bringing a variety of models under one unifying approach not only 
highlight the qualities of the quantum inverse scattering program but also allows us to better 
understand the relevant properties entering their Bethe Ansatz solution. This also motives us 
to look for further extensions which could shape our knowledge towards a possible classification 
of integrable models from an algebraic point of view. An interesting example seems to be the 
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vertex model, which we plan to investigate in a future work. 
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Appendix A : Boltzmann weights of the Shastry model 
We start this appendix by presenting the ten non-null Boltzmann weights of Shastry's 
i?-matrix (15). They are given by 

ai{X, /x) = {el'^('^)-''(^)la(A)a(//) + e-^''^''^-''^^H{X)b{ii)} a^iX, i^) (A.l) 

a2(A, /i) = {e-[''('')-'^(^)]a(A)a(/i) + e^''^^^-''^^H{X)b{ij)} a^iX, /x) (A.2) 

e['^('^)+^Wla(A)6(^) + e-[^(^)+^W]6(A)a(Ai) / cosh[/i(/i) - h{X)] \ , , . 

""'^^'^^ - a(A)6(A) + a(^)6(^) j cosh^^) + /.(A)] / "^^^'^^ ^^'^^ 

. e-['^('^)+'^Wla(A)6(^) + e[^(^)+^Wl6(A)a(/x) rcosh(/i(/x)-/i(A))l 

a(A)6(A) + a(^)%) \ cosh(M/.) + MA)) ' ^ ^ ^ 

[ e[MM)+MA)]a(A)fe(^) - e-tM^)+^Wl6(A)a(/x) ] ,,,, , ^,,,,, cosh[M/x) - MA)] . 

"^^^'^^ = I a(A)6(A) + a(^)6(^) / ^^^"^ ^^^^coshlM/.) + MA)]"^^^'^^ 

(A.5) 

[ -e-[M.)+MA)]^(A)fc(;,) +e[MMHMA)]fc(A)a(;,) | cosh[hi^^) - h{X)] 

"^^^'^^ = 1 a(A)6(A)+a(/.)6(/.) j^^ ^^^'^ ^^^KosHK^^) + h{X)f'^^^ 

(A.6) 

q;8(A, ii) = {e['^('^)-'^(^)la(A)6(//) - e-['^('*)-'^W]i,(A)a(//)} q;5(A, /i) (A.7) 

a9(A, /x) = {-e-['^('''-'^(^)]a(A)6(/i) + e['^('^)-''(^)l6(A)a(/i)} a5(A, /x) (A.8) 
X 62(^)_^2(A) r cosh[/i(/i) - /i(A)] 1 , 
"^'^^^'^^ = a(A)KA) + «(/.)&(/.) 1 cosh[M/.) + MA)] J ^^"^^ 
where the weight as (A, //) has been used as a normalization. We recall that functions a (A) and 

6(A) satisfy the free-fermion condition a^(A) + 6^(A) = 1, and in this paper we shall use the 
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parametrization is a(A) = cos(A) and b{X) = sin(A). There are certain useful identities satisfied 
by these weights we have used to simphfy commutation rules and the multi-particle problem. 
These relations are given by 



as^X, fi) = ai{X, fi) + aQ^X, fi) a^^X, fi) + a^i^X, fi) = a2{X, fi) (A. 10) 

tt2(A, /i)ai(A, /i) — ag{X, ix)as{X, jj) = a^^X, /x)a3(A, /i) — Q;io(A, /i) = (yl{X, jj) (A. 11) 
q;2(A, /i)a3(A, jj) + a4(A, /i)ai(A, /i) = 2a^{X, jj) (A. 12) 

Appendix B : Extra commutation rules 

This appendix is devoted to complement the commutation relations presented in the main 
text. For instance, there are some additional commutation rules which are important for the 
complete solution of the two-particle state problem. These are relations between the fields 
5(A), 5* (A), C{X) and C*{X) given by 

Ca{X)B,{^^) = -^^^B,{fi)Ca{X) + ^^^^[B{^l)A,,{X) - B{X)AM] (B-l) 
B:{X)B,{f,) = -^^|M5,(^)5*(A) + ^^^^[F{^,)AUX) - F{X)AM] (B.2) 

c:{x)B,if,) = ^^^Mf^)C;{X)-^^^B^^^^^^^ 

+^^4^^l^AUX)AM + ^ ''''^!''^h ab[F{^i)CiX) - B{fi)D{X)] 

(B.3) 

In particular, the commutation rule (5.3) is of considerable importance in the proof that 
the eigenvectors constructed in section 4 are highest weights states of the SU{2) rj- pairing 
symmetry (see section 5.2). In order to understand the role of the creation field B*{X) it 
is indispensable to derive its commutations relations with the other relevant fields. Between 
-B*(A) and the diagonal operators we have 

A{X)®B*{^^) = -i^lpfr(^,,X).[B*ifx)®A{X)]+t^^^^ 

a8(/i. A) aslAi, A) 

(B.4) 
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D{X)B- (rt = MD(X) - ,^^B'(X)D{i.) (B.5) 

as{X,iJ.) a8[X,iJ.) 

B(X)B-W = -,^2(jhRB^^^B(X)+'^F(u)C(X) 

'''^^'^'^F{X)C{fi)-t^4^[B{X)^A{^i)].e (B.6) 



aiiljL, X) a7(iji, X) 

while with itself and with the scalar operator F{X) we have 

S*(A)®5*(/x) = ^4^f(/x,A).[5*(/x)®5*(A)]+z^i^^{F(A)D (B.7) 
Qi2(A*, A) a7[ii,X) 

F(A)5*(/.) = ^^p^Ml^)B*iX) - z^^5*(/.)F(A) (B.8) 
a2[ii,X) q;2(/^, A) 

— * — * 

Lastly, the commutation rules with the annihilation fields C(A) and C*(A) are 

C:{X)B*M = -^4^i?,*(/.)C:(A) - t^^[D{f,)AUX) - D{X)AM] (B.IO) 
ttslA*, Aj a8[^,X) 

q;7(/x, A) arly", A) a?!/^, A) 

a7{ii,X) a7{ii,X) 

(B.ll) 

The best way of seeing that these later commutations relations are connected to those for 

— * 

the field B{X) is to read the equations in terms of their components. For instance, we note 
that commutation rule (B.2) is self-dual under the "dual" transformation described in section 
3. Several other relations have similar property as well. 

We close this appendix by presenting the expressions for the fundamental commutation rules 
when we solve the standard Yang-Baxter algebra (3). These relations lack the presence of the 
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imaginary factors "i" and certain extra signs when compared to their graded counterparts. 
Below we hst the most important relations for the creation fields B{X) and F{X) 



^\b{i,) ® A{X)].fTw{X, 1^) + ^4M^(A) ® Ail,) 



aio(A, /x) 
a7{X,n) 



B{X)B{ii) 



ag{X,ny 

[B*{X)B{i,) - ^4Mf(A)C(/.) + ^^%I^F{i,)C{X)]®(tw 



ag(X,iiy 



ag(X,iJ,y 



agdi^X)' 



ag{ijL,Xy 



(B.12) 
(B.13) 



D{X)B{t,) = 



'''^^^^'^-mDiX) - ^^F{u)C*{X) 



a7{X, fi) a-j{X, jj) 



(B.14) 



Aab{X)F{ll) 



ag{X, Id) as{X, ijl) 

"5(A,/i) 



]F(/.)i„,(A) - 



n?(A./0 



ag{X, ljL)a8{X, ijl) 



F{X)AM 



— 7T—Am^B*{t,)],^+^^^^^[^^ 

ag{X,ll) 0(8{X,ll) 



(B.15) 



5(A)F(;.) = -^^Fii.)BiX) + ^^FiX)Bii.) + '^^^ (B.16) 



D(A)F(/x) 



a7{ii,x) ^^{n.xy 

"^^^'^V(^)Z^(A) + ^#4^(A)^(/^) 



q;7(A,/x) Qi7(A,At) 

— * 

where ^tvk and r(A, are given by 



a'j{ii, X) 

aio(A,;u) 
Qi7(A,/x) 



W.{5*(A)®5*(/i)} 



(B.17) 



eriy= (0 110) fTi^(A,//) 



(B.18) 



/I 0\ 

a(A,/x) -b{X,n) 

-b{X,fi) a(A,/i) 

Vo ly 

Furthermore, the relations closing the commutation rules between the creation operators 
B{X) and F(A) are 

B{X) ® B{ix) = ^4M[^(/^) ® B{X)].fTwiX,f^) - ^4^{F{X)B{id) - F{ix)B{X)}iTw 



q;2(A,//) 



q;7(A,//) 



(B.19) 
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F{X)B{i,) 



[F(A),F(/.)] = 



Qi2(A,/i) 
«5(A,/i) 



fi(^)F(A) 



tt2(A,/i) 

a9(A, /i) 



F(/x)E(A) 



(B.20) 
(B.21) 

(B.22) 



Q;2(A,yu) ' ' a2(A,^) 

Appendix C : The two-peirticle state 

In this appendix we provide details about the technical points entering the analysis of the 

two-particle eigenvalue problem. We begin the discussion by first considering the wanted terms. 

We recall that the amphtudes proportional to the first part of the two-particle eigenstate are 

easily estimated as a product of the first right-hand side terms of the commutation rules (34- 

36). For the second part, however, there are more contributions since the action of diagonal 

operators on the first part -B(Ai) (8) -B(A2).^ |0) produce at least one extra term proportional to 

— * — * 

the second part F(Ai)^.J^|0) as well. It turns out, however, that these contributions miracu- 
lously factorize in the same product forms we have obtained for the first part of the eigenstate. 
This happens thanks to remarkable identities between the Boltzmann weights we begin listing 
below. For the field B{X) there are two contributions and they factorize as 

a2{y,x) a2{y,x) a5{z,x) aio{y,x) a7{y,z) a2{y,x) a2{z,x) 



a7{y,x) a9{y,x) ag{z,x) a7{y,x) aio{y,z) ag{y,x) ag{z,x) 

Analogously, for the field -D(A) we have 

a2{x,y) _ a-,{x,y) aio{x,z) aijy^z) ^ as{x-,y) asix^z) 

a7{x,y) a7{x,y) a7{x,z) aio{y,z) a7{x,y) a7{x, z) 



(C.l) 



(C.2) 



For the diagonal field ^ Aaa(A) we have three contributions, where two of them are gen- 

a=l 

erated by the first part of the eigenstate. The identity that brings these terms together and 
also gives rise to the auxiliary eigenvalue function is 



1 + 



al{x,y) 



as{x,y)ag{x,y) 



ai{x, y) aio{x,z) a5{x,y) a7{y,z) _ 

—( 7 -( T—( 7 1 7 1 + « y)\ 

ag[x,y) a7{x,z) a8{x,y) aio{y, z) 
^ aio{x,y) a2ix,y) a^jx, z) a7{y,z) 

a7{x,y) ag{x,y) ag{x,z) aio{y,z) 

ai{x, y) ai{x, z 



ag{x,y) ag{x,z) 



b{x, y) + b{x, z) — a{x, y)a{x, z) 



(C.3) 
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— * — * 

Next we turn to the analysis of the unwanted terms proportional to B{X) B{Xj) and 

— ♦ — ♦ A — * 

(8) {B*{X) (8) I)] (8) B{Xj). The terms with Xj — X2 are straightforwardly read from the 
commutation rules (34-36) because only single contributions occur for each diagonal field. 
However, for Xj = Ai, the situation is more complicated because it involves many different 
contributions whose origin is due to the fact that the rapidity Ai is wrongly ordered when 
compared with A2. Nevertheless, one expects that there should be a better way of recasting 
these terms since the Bethe Ansatz equations are usually independent of indices relabeling. 
Indeed, it turns out that these many contributions can be compactly written by introducing 
the "ordering" factor Oj^\X,Xj; {A^}). As before, in order to factorize these contributions to 
a single term, we had to use extra identities between the Boltzmann weights. For example, for 
the field B{X) they are 

—f ^—f ^a{y,x) — -— - = — -— ra{y,z) (C.4) 

aQ{y,x) ag{z,x) a7{y,x) a7{y,z) aQ{y, x) ag{z,y) ag{z, x) ag{y, z) 

and 

aw{y,x) aio{y,z) ai{y, x) a^iz, x)- a^iz, x) ai{y, z)- 

a7{y,x) a7{y,z) ag{y, x) ag{z, x) ag{z, x) ag{y, z) 

where the left-hand side of the above equations represent the contributions coming from the 

"brute force" calculations while the right-hand side exhibits the "ordering" factor explicitly. 

2 

Similar simplifications can be carried out for the fields ^^ao(A) and -D(A), but we skip 

a=l 

further details since there is a much simpler way to understand the origin of such "ordering" 
factor. As it has been explained in section 4, this factor can be easily derived with the help of 
the exchange property (69). Anyhow, the coincidence between the "brute-force" computations 
and the symmetrization results gives us confidence to go ahead using the symmetrization 
procedure for multi-particle states. 

Finally, we show how the third type of unwanted terms generated by the diagonal field 

2 

^Aaa(A) can be further simplified. First it is convenient to rewrite the term proportional 

a=l 

to [a;i(Ai)a;2(A2)]^F(A)^.^ in a way that the auxiliary eigenvalue function appears explicitly. 
For this purpose we use the second identity (66), and rewrite the contribution to the above 
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mentioned unwanted term as 



K(Ai)w2(A2)]^A(^)(A = A2, {Xi})H-i{X, Ai, A2)[6(Ai, A2) - a(Ai, X2)]F{X)i.f (C.6) 

Next we take advantage of the symmetrization property of the two-particle eigenstate and 
evaluate the contribution proportional to [uji{X2)uj2{Xi)] F{X)^.J^ as follows. The idea is to 
begin with the right-hand side of equation (69), which remarkably gives us precisely the extra 
r-matrix necessary to produce the auxiliary eigenvalue at A = Ai. Obviously, the amplitude 
contributing to this term is proportional to function obviously H^IX, A2, Ai) multiplied by the 
extra factor coming from the exchange relation. Putting these information together 

we are able to rewrite the second contribution as 

- K(A2)^/;2(Ai)]^A(i)(A = Ai, {Xi})H,{X, A2, X^ f'[^'' ^'I f{X)It (C.7) 

These manipulations make the cancellation of the third type of unwanted terms more 
transparent, since it allows us to use the Bethe Ansatz equations in a more direct way. Indeed, 
using the Bethe Ansatz equations (65) in the terms (C.6) and (C.7) and adding them to 
those coming from the fields B{X) and D{X), we find that the unwanted terms proportional to 
F(A){.JF are cancelled out thanks to the following identity 

Hi{x,y,z) + H2{x,y,z) = H^{x,y, z)^{y, z) - a{y,z)] - H^ix, z,y)^^^j^ (C.8) 

a2{y,z) 

This gives us another opportunity to verify the symmetrization scheme. Comparing (C.8) 
and (67) we conclude that the following identity 

(X\ (y^ z^ — 
H3ix,z,y) ' = H4{x,y,z)[a{z,y) - b{z,y)] (C.9) 
Oi2[y,z) 

is indeed satisfied. 

We remark that the above technicalities are of enormous help when we consider generaliza- 
tion to multi-particle states. In next appendix we shall discuss this fact for the three-particle 
state. 

Appendix D : The three-particle state 
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We shall start this appendix showing how the permutation symmetry Ai <-> A2 is imple- 
mented for the three-particle state. As before, our strategy consists in reordering the rapidities 
Ai and A2 with the help of the commutation rule (25). This allows us to write the Ansatz (76) 
as 



$3(Ai, A2, A3) 



a2(Ai, A2) 
.a;io(Ai, A2) 

1 

"7(^1, A2) 



tt7(Al, A2 

[FiX2)B{X^)i^ B{Xs)] + i^^^2lM[B{X^) ® iF{X2)B{Xs)] 



o;7(A2, A3) 



+[e ® F(Ai)i?(A3)i?(A2)]#(Ai, A2, A3) + [e ® F(Ai)i?(A2)i?(A3)]^r (Ai, A2, A 

Next we use the commutation rule (35) to simplify the second and the third parts of the 
above equation, earring the scalar field B{Xj) {j = 1,2) through the creation operator BlX^). 
This procedure not only helps us to eliminate the fifth term of equation (D.l) but also prompts 
the appearance of a desirable term proportional to (g) F(A2)-B(A3)S(Ai)]. Now, imposing the 
exchange property (78) for the rapidities Ai and A2 we find the following necessary condition 

tt2(A3, Ai) 



(3), 



(D.l) 



-(3)f^ X X N «7(A2, Ai) 



(D.2) 



aio(A2,Ai) q;9(A3, Ai) 
This relation together with the previous restrictions found in section 4, cf. equations (73- 
74), are able to determine unambiguously the constraints for the three-particle state. The 
next step is to show the consistency of the whole procedure, i.e. that the equality between the 
remaining terms are indeed satisfied. By using the commutation rules (41-42) we derive two 
consistency conditions, given by 

aio(A2, A3) a5(Ai, A2) tti(Ai, A2) aio(Ai, A3) 



[B{X2)^C]F{X^)B{Xs) 
[i0B{X2)]F{X,)B{Xs) 



a7(A2,A3) a2(Ai,A2) a2(Ai,A2) «7(Ai, A3) 



■n2(Ai, A2) 



aio(Ai, A2) q;5(A3, A2) a8(Ai, A2) , «8(Ai, A2) .(3).. . . s 

H 7^ — r^92 (Ai, A2, A3) 



"7(Ai,A2) a9(A3, A2) a2(Ai, A2) a2(Ai,A2)' 



(D.3) 



and 



F{X2m®B{X^)]B{Xs) 



Q;io(Ai, A2) q;5(A3, Ai) a5(Ai,A2 



+ 



#(Ai,A2,A3) 



q;7(Ai,A2) q;9(A3, Ai) q;2(Ai,A2) 
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aio(^i, ^2) "5(^3, ^2) «5(Ai, A2) ai(Ai, X2) .(3) 



Qi7(Ai,A2) q;9(A3, A2) a2(Ai, A2) a2(Ai,A2; 



92 (A2, Ai, A3).fi2(Ai, A2) 



= -F(A.)[B(Ai) »^1i?(A3) °'°'''-'°' If'- (D.4) 

«7(A2,A3) Q;2(Ai,A2) 

In order to disentangle the above expressions we need the help of certain useful identities 
between the "exclusion" vector and the auxiliary r-matrix. More precisely, they are 

[e ® B{y)]h2{X, f^) = [a(A, fx) - b{X, B{y)] (D.5) 

[e®5(|/)]f23(A,/i) = [i0B{y)] + b{X,f,)[B{y)®i] (D.6) 

[B{y) e>23(A, //) = [a(A, fx) - b{X, //)] [B{y) ® (D.7) 

[^(|/)®e>i2(A,/x) = [5(|/)®cl + 6(A,/.)[e®S(|/)] (D.8) 

Inserting the identities (D.5-D.8) into equations (D.3-D.4) we end up with four identities 
among the Boltzmann weights which have been verified by using Mathematical'^ . With this 
we complete the symmetrization analysis for the three-particle state. 

We now turn to the analysis of the eigenvalue problem for the thrce-particlc state. Let 
us begin by investigating the action of the scalar field B{X) on the state (76). The first step 
consists to carry the field B{X) through the creation fields -B(Ai) and -F(Ai) by using the 
commutation rules (35) and (41-42). Afterwards, we use directly the known results for the 
two-particle state, cf. (57), in order to turn one more time the scalar fields B{X) and -B(Ai) 
over the two-particle state |$2(A2, A3)). As a third step, we need to reorder creation fields such 
as -B(Ai) and B{X) with the help of commutation rule (25) as well as keep on carrying the 
scalar field B{X) until it reaches the vacuum. After this long but straightforward computations 
we find the following result 

S(A)|$3(Ai,A2,A3)) = ki(A)]^ n ^ ""'i^' ^\ B{X,) ^ $2{X2, Xz).^ |0) 

j=i "i(Aj,A) 

+ki(A)]^i^4T^^[^' ® F{X,)B{Xs)B{X2)g['\X^, A2, Xs).f |0) 

a9(A3, A) 

+K(A)]^i^4T^^[^' -^1' -^3]^ ® F(Ai)B(A2)i?(A3)#(Ai, A2, A3)./ 10) 
a9(A2, Aj 

-|-unwanted terms (D.9) 
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where function L[x,y,z] is precisely the left-hand side of identity (C.l) we have worked out 
for the two-particle state. This allows us to factorize the amphtudes for the second and the 
third terms of the above equation, and as result all the three wanted terms have a common 
amplitude as it should be. To what concerns the unwanted terms, our computation shows that 
they can be gathered in two basic families. More specificly, they are proportional to 

K(A,)]^^(A)0$2(A,,Afe) (D.IO) 

and 

K(A,)«;i(AO] F(A)$i(Afc) (D.ll) 

The first term in the family (D.IO), say Xj = Ai, A; = A2 and A^ = A3, is originated from 
the first part of the three-particle state when we turn the scalar field B{X) through B{Xi). 
Keeping the second term of the commutation rule (35), and by using the two-particle results 
(57) to carry -B(Ai) through |$2(A2, A3)), we find that its amplitude is 

Q;9(Ai, A) Q;9(Ajfc, Ai) 
We estimate the amplitudes of the remaining terms in the family (D.IO) by taking into 
account the exchange property (78), in much the same way we did for the two-particle state. 
This means that the amplitudes are going to be multiplied by the first "ordering" factors 
Of\Xj; {Afe}), and three possible unwanted terms j — 1, 2, 3 can be compactly written as 

-K(A,)]^z^4r4 n '^^^^AX)0^,{X^, . . . , A„ . . . , A3) xO«(A,; {A,}).^ |0) (D.13) 

The contributions to the second family of unwanted terms come from all the pieces com- 
posing the three-particle state. It turns out that for k = 2,3 their amplitudes can be computed 
in a very similar way we did for the second and third parts of the wanted terms, respectively. 
The main difference is that now we have to keep track of terms proportional to F{X) rather 
than -F(Ai). We find that the amphtudes for these unwanted terms are 

K(A>,(A,)]^i/i(A,A„A,) n ^?7r44''?7I^4 [nA)f®$i(A2,...,A,,...,A3)] 

xdg.)(Ai,A,;{Afe}).J?|0) (D.14) 
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where 0^j{Xi, Xj] {Xk}),j — 1, 2, is the second type of "ordering" factor which has been already 
defined in the main text, see equation (88). It should be emphasized that we have derived the 
above factor from a "brute force " analysis, and similar to what happened to the two-particle 
state, this gives us the clue to proceed in order to better estimate the remaining unwanted 
terms appearing in this family. We easily recognize that this factor is related to the operation 
of bringing two rapidities in the first two positions of the eigenvector. Keeping this in mind, 
we see that all the contributions to the second family of unwanted terms can be written by 

[w^{Xi)w^{Xj)]'^H,{X, A,, A,) n ^^^^^^^ ^^(^i' ■ ■ ■ , A,, . . . , A„ . . . , A3)] 

t=i "9(Afe, A;) ag{Xk,Xj) 

x6;f(A,,A,;{Afe}).:?|0) (D.15) 

Collecting the expressions (D.9),(D.13) and (D.15) we find that the action of the scalar field 

B{X) on the three-particle state is described by the formula (82) with n — ?>. Similar reasoning 

2 

can be repeated for the fields D{X) and ^ Aaa{X)^ and only when we are estimating the third 

a=l 

type of unwanted terms new technicalities emerge. In what follows we present the details of 
these computations in the simplest case, i.e. the situation where no "ordering" factors are 
needed. Generalization for the remaining terms is along the lines of formula (D.15). For the 
field -D(A) we find that such amplitude is 

//2(A, Ai, A2)F(A)6cA6i(Ai)ic6.(A2)Sb3(A3):^'^'^'^ |0) (D.16) 

Now, carrying the operators ^^^^(Ai) and Acb^iX^i) through B},^{X^) with the help of commu- 
tation rule (34) we find 

K(Ai)^2(A2)]^//2(A,Ai,A2) n ^^\ i^cF{X)B,{X,y^^-{X,, \^)rt{>^i. X,)r'-'^'^ |0) 

ZQ!9(Afe, A3) ^ 

(D.17) 

which is further simplified by using the following identity 

T(i)(A = A2, {A,})SS|r^'nA = Ai, {A,})--^3 = ^M3(^^^ A3)f;r(Ai, A3) (D.18) 
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Inserting (D.18) into (D.17) we finally obtain 

[w2{Xi)w,{\,)]'^H,{X, Ai, \2)A^'\X = Ai, {A,})A«(A = A^, {A,}) n ^li^^hM_F{\)^® ^(Aa) 

^J-^ tag[Ak, As) 

(D.19) 

2 

For the field ^ Aaa{X) we find that one of the contributions is 

a=l 

i^3(A,Ai,A2)U^(A)5(Ai)i„fc,(A2)56,(A3)^'^'^'M0) (D.20) 
and when we carry -B(Ai) and Aab2{Xi) through ^^^(As) we have 

[^i(Ai)^2(A2)]^i^3(A,Ai,A2) n n'' \'\ uF{X)f',f-^{X,, X,)B,{Xs)J''^'-'^ |0) (D.21) 
Next using the following identity 

UAt'iX^, X3)B,{X,)T'-^'-'^ = e,5f^f (Ai, A2)T(^)(A = A^, {Xi}t'pt B,{X,)T'-^'-'^ (D.22) 
we finally find 

K(Ai)«;2(A2)]^i/3(A, Ai, A2)[a(Ai, A2)-6(Ai, A2)]A«(A = A2, {A,}) fl ^^^j^^F{X)i®B{X;) 

k=i «a9(Afc, A3) 

(D.23) 



Lastly, the second contribution coming from the field ^ ^^^(A) is estimated by using the 

a=l 

same trick explained in the previous appendix for the two-particle state. We further remark 
that the technical points explained in appendices C and D are valid for many other models 



such as the Bariev XY chain and those solved in ref. . 

Appendix E : The Bariev model 

We start this appendix by presenting the Boltzmann weights pj{X,fi)j = 1, . . . , 15 found 
by Zhou p7|. Normalizing them by pi{X,fi) we have 



P3(A,p) = 1^^12x1^ ' P4(A,/i) = p2(A,/i)p2(A//i,/i//i) (E.2) 
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1 



1 + h^XfjL 



+ 



(l + Xpi) 



;i + AAt)(l + /i2A/x) 

P3(A,^)(A - /iV) 



h 



, ps{\,p) = P2i\/h,fi/h) 



Pio(A,/x) 



+ 



1 + A^ + /i^ - A/x 



l + /i2A/x (1 + A/x)(l + /i2A/x) 



, pii(A,/x) = p3(A//i,/i//i) 



Pl2(A,//) = -P5(/^, A) , Pl3(A,//) 



Pl2{\p) 
h 



(E.3) 

(E.4) 

(E.5) 

(E.6) 

(E.7) 
(E.8) 



Pi4(A,/x) = p9(/x, A) , pi5(A,/x) = p3(A,/i)pii(A,/x) 

We remark that a rescaling A ^ and /i ^ brings these weights in a more symmetrical 
form and this is useful, for example, to check the Yang-Baxter equation (6). Besides that there 
are some few identities between those weights that were useful in the calculations of section 6. 
They are 

Pi5(A, p)[p9{X, p)+pi(X, /x)] = P5(A, /x)p6(A, /x), P6(A, p)pi(X, p)+P5(X, /x)pi5(A, /x) = p6(A, /x)p7(A, /x) 

(E.9) 

Pi2(A,/x)[p9(A,/x)+Pi(A,/x)] = P5(A,/x)p4(A,/i), p5(A,/x)pi(A,/i)+pi5(A,/i)p6(A,/x) = P5(A, /i)pio(A, /x) 

(E.IO) 

Next we present the commutation relations we found by solving the standard Yang-Baxter 
algebra (3). Between the diagonal fields and the creation operators we have 

i(A)0^(/x) - ^^[B(/x)®i(A)].f(^)(A,/x)-^^B(A)®i(/x) + 



P^{\p) 
Pb{X,n) 
P9{X,n) 



PsiX,^)' 



PsiKp)' 



Pz{\p)Pb{X,p) 



psipAY 



psipAY 



(E.ll) 
(E.12) 



D{X)B{^^) = ^4M^if^)DW + ^%4nP)C*W 



P9{X,p) 
Pi{\p) 



P9{X,p) 



FiX)C*{i,)-P^^e''^.{B*iX)®A{i,)} 

P9{X,P) P9{X,P) 



(E.13) 
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P3(A,yU) Pll(A,/i) p3(A,/i)pii(A,//) 

[S(A) + f^i^[B*{X) ® (E.14) 



P9(/i, A) P9(/i, A) P9(/^, A) 



i^(A)F(;.) = ^4^F(;.)i^(A) - ^4M,.(A)i^(;.) - ^4^e^-).{i?*(A) ^ (E.16) 

P9[A,H) P9[\P) P9[\P) 

Between the creation fields we have 



B{X)^B{fi) = [S(/x)(g)S(A)].f(^)(A,/x) + 



P9{\P) P9{\P) 



(E.17) 

[F(A),F(/x)] = (E.18) 
F(A)5l;.) = '-^,Fi,)BiX) + '-^mF{X) (E.19) 

5(A)F(/.) = ^4^5(;.)F(A) + P^j^Fi^^)BiX) (E.20) 
The extra relations for the analysis of the two-particle state are 

Ca{X)B,{i,) = B,{ii)Ca{X) - P^!^[B{X)Aa,{ii) - B{ii)Aa,{X)] (E.21) 

PsiXP') P3(A,//) P3(A,//) 

C:{X)B,{^^) |0) = ^4r4 [^('")^(^) - ^-(^)^''''(/^)] |0) (E-23) 

P9(A,//) 
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